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Abstract. We present our image processing system for the reduction of op tical imaging data fr om multi-chip 
cameras. In the framework of the Garching Bonn Deep Survey (GaBoDS; ISchirmer et ai] |2003^ consisting of 
■ about 20 square degrees of high-quality data from WFI@MPG/ESO 2.2m, our group developed an imaging 

' pipeline for the homogeneous and efficient processing of this large data set. Having weak gravitational lensing as 

the main science driver, our algorithms are optimised to produce deep co-added mosaics from individual exposures 
obtained from empty field observations. However, the modular design of our pipeline allows an easy adaption to 
different scientific applications. Our system has already been ported to a large variety of optical instruments and 
its products have been used in various scientific contexts. In this paper we give a thorough description of the 
algorithms used and a careful evaluation of the accuracies reached. This concerns the removal of the instrumental 
signature, the astrometric alignment, photometric calibration and the characterisation of final co-added mosaics. 
In addition we give a more general overview on the image reduction process and comment on observing strategies 
where they have significant infiuence on the data quality. 
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1. Introduction 

During the last decades, optical Wide-Field Imaging has 
become one of the most important tools in observational 
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astronomy. The advances in this field are closely linked 
to the development of highly sensitive detectors based 
on charge coupled devices (CCDs). These detectors have 
grown in size from a few hundred pixels on a side to the 
currently largest arrays of about 4000 x 4000 pixels. The 
need for ever larger fields-of-view and deeper images, and 
the technical constraints in the manufacturing processes 
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Fig. 1. Shown is the CCD layout of WFI@MPG/ESO 2.2m (not to scale). The mosaic design of CCDs to obtain a larger field- 
of-view generally confronts the user with new issues compared to the treatment of single-chip cameras. The possible overlap of 
sources lying on different detectors on dithered exposures significantly complicates the astrometric and photometric calibration 
and the subsequent co-addition process. First, the large field-of-view leads to notable optical distortions from the field centre 
to the edges (see also Fig. IIH . These distortions have to be known accurately for each CCD to create a correctly stacked 
mosaic. Second, for homogeneous photometry we need to take into account different quantum efficiencies and hence different 
photometric zero points of individual chips. Finally, dithering leads to non-uniform noise properties in the co-added image (due 
to the gaps, sensitivity and gain differences; see also Fig. l24l l. The knowledge of the relative noise properties between individual 
pixels in the co-added mosaic is important in the subsequent source catalogue creation process (see Fig. 1251 . 



of larger CCDs led to the development of multi-chip cam- 
eras in the mid 90's. Hereafter, we refer to them simply 
as Wide Field Imagers, or WFIs. The production of such 
a mosaic detector array with well aligned CCDs in three 
dimensions is very difficult. The misalignments of indi- 
vidual chips, and the multi-chip nature by itself, are the 
main reasons that the data reduction process is signifi- 
cantly more complicated compared to the treatment of 
single-chip cameras. Although a multi-chip CCD camera 
seems to be a simple collection of N single-chip cameras 
(see Fig.nj at first glance, it is, in general, not possible to 
treat all images belonging to a certain detector indepen- 
dently during the complete reduction process. -'^ To obtain 
an accurately aligned co-added image on the whole WFI 
scale we have to properly account for the mosaic structure 
with its gaps, not perfectly aligned chips, and a large field- 
of-view that needs to be treated as a curved surface instead 
of a flat plane. We have to apply sophisticated techniques 

^ This is still possible for certain scientific objectives not 
needing any image co-addition (as variability studies in single 
frames for instance) or if data are obtained in a very compact 
dither pattern so that different chips never overlap. As we will 
show in the course of the paper, compact dither patterns have 
in general severe drawbacks on the final image quality. 



going considerably beyond the work that has to be done 
for a single-chip camera with a small field-of-view. In ad- 
dition, a homogeneous photometric calibration and the 
characterisation of the noise properties in a co-added mo- 
saic poses new challenges on calibration techniques. Aside 
from the technical issues related to the mosaic structure of 
new cameras, the increasing quality of new optical instru- 
ments with their high sensitivity and throughput demand 
the use of state-of-the-art algorithms to optimally exploit 
the obtained data scientifically. 

Up to now, no standard procedure to process and treat 
WFI data is established among observers. The main rea- 
son is that WFI capabilities have only been available to 
a broader community for about half a decade and we 
are still in the process of fully understanding the prop- 
erties of the new instruments and developing the neces- 
sary tools for the data handling, and also for analysing 
reduced images. Another issue is the high data flow and 
the demands on computer hardware connected already to 
small and medium-sized WFI projects. Only for about two 
years high-end, and for the occasional observer affordable, 
desktop PCs are equipped adequately. 

That these issues are non-trivial is illustrated by the 
relatively small number of publications based on WFI 
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projects compared to the amount of data acquired. For 
example, after the WFI at the 2.2m MPG/ESO telescope 
(hereafter WFI@2.2m) became operational in January 
1999, only 35 refereed papers based on its data appeared 
until the end of 2002. Meanwhile, the rate has risen^, 
but is still significantly behind those of other ESO instru- 
ments which are under similar pressure by proposers as 
WFI@2.2m. 

In this paper we present the methods and algorithms 
we use to process WFI data. The presentation is organised 
as follows: In Sect.Qlwe introduce the basic concepts of the 
GaBoDS pipeline and the philosophy behind the choices 
we made. We discuss the advantages and the disadvan- 
tages that arise thereof. The pre-reduction process (i.e. 
the removal of the instrumental signature from the data), 
which is in principle identical for single- and multi-chip 
cameras, is described in Sect.O Details of the astrometric 
and photometric calibration of the data are presented in 
Sect.0] An explanation of our adopted scheme to deal with 
the inhomogeneous noise properties in co-added mosaic 
data is given in Sect.El followed by the image co- addition 
methods in Sect.|Hl We perform quality control checks on 
CO- added data in Sect. |7| and draw our conclusions in Sect. 

El 

No astronomical pipeline will produce the best possi- 
ble result with data obtained in arbitrary conditions and 
strategies. Hence, besides a pure description of algorithms, 
our presentation also contains guidelines how WFI data 
should be obtained to achieve good results. 

We assume familiarity with data processing of opti- 
cal imaging data. In this publication we focus on describ- 
ing the algorithms used for all necessary processing steps. 
Where the reduction differs significantly from well estab- 
lished algorithms for single-chip cameras (this mostly con- 
cerns the astrometric alignment and the photometric cal- 
ibration) we give a scientific evaluation and a thorough 
estimation of the accuracies reached. 

We note that we mainly work on data from WFI@2.2m 
and most of the examples and figures in this paper refer 
to data from this instrument. The reader has to be aware 
that the quoted results, the accuracies and the overall us- 
ability of our algorithms can vary significantly when being 
applied to data sets from other cameras. At critical points 
we come back to this issue in the text. 

2. Pipeline characteristics 

2.1. Scientific motivation 

While all optical data need the same treatment to remove 
the instrumental signature (bias correction, flat-fielding, 
fringe correction etc; see Sect. EJ, the subsequent treat- 
ment of the images strongly depends on the scientific ob- 
jectives and on the kind of data at hand. Our primary 
scientific interests are weak gravitationa l lensing studies 
(see e.g. iBartelmann fc Schneidedf200ll for an extended 

^ The ESO Telescope Bibliography turns up a total of 107 
refereed papers until December 2004. 



review) in deep blank-field surveys. These studies mainly 
depend on shape measurements of faint galaxies. To en- 
sure that the light distribution is deformed as little as 
possible by the Earth's atmosphere and the optics of the 
telescope, weak lensing data are typically obtained under 
superb seeing conditions at telescopes with state-of-the- 
art optical equipment. For a proper treatment of those 
data, the main requirements on a data processing pipeline 
are the following: 

— We need to align very accurately galaxy images of 
subsequent exposures that are finally stacked. This in- 
volves an accurate mapping of possible distortions (see 
Fig.im. 

— We need the highest possible resolution, and hence a 
co-addition on a sub-pixel basis (see Sect.EJ. Together 
with the previous step, this is crucial for an accurate 
measurement of object shapes in the subsequent lens- 
ing analysis. 

— We need to accurately map the noise properties in our 
finally co-added images (see Fig. I25|l . This knowledge 
enables us to use as many faint galaxies as possible and 
to push the object detection limit. This also requires 
that the sky-background in the co-added image is as 
fiat as possible. 

The algorithms we use are chosen to go from raw im- 
ages to a final co-added mosaic with the objectives de- 
scribed above. The responsibility of the pipeline ends with 
the co-addition step. We are fully aware that the cho- 
sen procedures and algorithms may not be optimal for 
projects having different scientific objectives such as ac- 
curate photometry, the investigation of crowded fields or 
the study of large scale, low surface brightness objects 
for instance. However, the modular design of our pipeline 
(see below) allows an easy exchange of algorithms or the 
integration of new methods necessary for different appli- 
cations. The main characteristics of our current pipeline 
are summarised in the following: 

— Ability to process exposures from a multi-chip camera 
with N CCDs. To date we have successfully pro- 
cessed data from: WFI@2.2m, CFH12K@CFHT, 
F0RS1/2@VLT, WFI@AAO, MOSAIC- 
I /II@CTIO /KPNO , SUPRIMEC AM@Subaru, 
WFC@INT, WFC@WHT, and several single-chip 
cameras (e.g. BUSCA@2.2m Calar Alto). 

— Ability to handle mosaic data taken with arbitrary 
dither patterns. 

— Precise image alignments without prior knowledge 
about distortions. The astrometric calibration is per- 
formed with the data itself. 

— Absolute astrometric calibration of the final images to 
sky coordinates. 

— An appropriate co-addition of data obtained under dif- 
ferent photometric conditions. 

— Image defects on all exposures are identified and 
marked before the co-addition process. 
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- Creation of weight maps taking into account differ- 
ent noise/gain properties and image defects for the co- 
addition. 

- Accurate co-addition on sub-pixel level/ability to 
rescale data to an "arbitrary" scale (ability to easily 
combine data from different cameras/telescopes). 

We are currently extending our algorithms to near-IR 
cameras which will be described in a forthcoming pub- 
lication (Schirmer et al. in prep.). 

2.2. Implementation details 

As mentioned above, many tools for WFI data process- 
ing are currently under active development and different 
groups have already released excellent software packages 
for specific tasks. Hence, we built our own efforts on al- 
ready publicly available software modules wherever possi- 
ble. Many of the algorithms used are very similar to those 
developed for the EIS (ESQ Imaging Survey) Wide survey 
in 1997-1999 fsee lNonino et al.lll999|) . The main pillars of 
our pipeline are the following software modules: 

- The LDAC software^: The LDAC (Leiden Data 
Analysis Centre) software was the backbone of the first 
EIS pipeline. It provides a binary catalogue format (in 
the form of system-independent binary FITS tables) 
including a large amount of tools for their handling. 
Moreover, this module contains software for the astro- 
metric and photometric calibration of mosaic data. 

- The EIS Drizzl e^: A specific version o f the IRAF 
package drizzle (|Fruchter fc Hoo3l2002l) was devel- 
oped for EIS. It directly uses the astrometric calibra- 
tion provided by the LDAC astrom part and performs 
a weighted linear co-addition of the imaging data. 

- TERAPIX software^: (Bertin et al. 2002) 
SExtractor is used to obtain object catalogues 
for the astrometric calibration. Moreover it produces 
a cosmic ray mask in connection with EYE^ in addition 
to smoothed and sky-subtracted images at different 
parts of the pipeline. SWarp, the TERAPIX software 
module for resampling and co- adding FITS images, is 
used alternatively to EIS Drizzle for the final image 
co-additiom 

- FLIPS'^: l|Magnier fc Cuillandr(j2004j) FLIPS is one of 
the modules for data pre-reduction. It is optimised to 
perform operations on large format CCDs with min- 
imal memory requirements at the cost of I/O perfor- 
mance. 

^ available at |ftp: / /ftp . strw.leidenuniv.nl/pub /Idac / software / f 
available via ESQ 's SCISOFT CD (see ' 

|http:/ /www.eso.org/scisoft / \ 

available at http://terapix.iap.fr/soft/ 
® EYE (Enhance Your Extraction) allows the user to generate 
image filters for the detection of small-scale features in astro- 
nomical images by machine learning (neural networks). These 
filters are loaded into SExtractor that performs the actual de- 
tection. We use such filters to detect cosmic rays in our images. 
see ,http://www. cfht.hawaii.edu/~jcc/Flips/flips. html, 



Eclipse and qfits tools": llDevillardI l200lh We 

use several stand-alone FITS header tools from the 
Eclipse package to update/query our image header. 
Moreover, tools based on the qfits library offer an 
alternative to FLIPS for image pre-reduction. 

— Astrometrix^: Astrometrix (developed at 
TERAPIX) is another module for obtaining as- 
trometric calibration. 

- IMC AT utilities^": From the IMCAT tools we exten- 
sively use the image calculator ic. 

These tools have been adapted for our purposes if neces- 
sary and wrapped by UNIX/bash shell scripts to form 
our pipeline. Our main effort is to provide the neces- 
sary interfaces for the communication between the indi- 
vidual modules and to add instrument and science spe- 
cific modules for our purposes (as a thorough quantifica- 
tion of PSF properties for instance). With our implemen- 
tation approach we can make use of already well-tested 
and maintained software packages. The description above 
also shows that we can easily exchange modules as soon 
as new algorithms or better implementations for a cer- 
tain task become available. To further enhance the mod- 
ularity of the pipeline, we split up the whole reduction 
process into small tasks (a complete reduction process for 
WFI@2.2m data typically involves the call to 20-30 dif- 
ferent scripts but superscripts collecting several tasks can 
easily be written). This ensures a very high flexibility of 
the system and that potential users can easily adapt it to 
their needs or to specific instruments. 

Nearly the whole system is implemented in ANSI C 
and UNIX bash shell scripts (EIS Drizzle is written in 
FORTRAN?? and embedded into IRAF, Astrometrix is 
developed under Perl-fPDL and parts of our photomet- 
ric calibration module are implemented in Python) . Thus 
we have full control over the source codes and can port 
the pipeline to different UNIX flavours. Table lists the 
systems on which we successfully used it so far. 

The main disadvantage of building up a pipeline from 
many different software modules instead of developing 
a homogeneous system from scratch is that it becomes 
very difficult to automatically control the data flow or 
to construct a sophisticated error handling and data in- 
tegrity checking system. Furthermore, our pipeline so far 
offers only very limited possibilities to store the history 
of the reduction process or to administrate raw and pro- 
cessed image products. Also, formal speed estimates for 
the throughput of our processing system are low compared 
to homogeneous systems (see Table EJ- Thus, the usabil- 
ity for large, long-term projects of the system is limited 
at this stage. 

With the compactness, the flexibility and the usabil- 
ity of our system for the occasional user^^ we regard 

* available at http://www.eso.org/projects/aot/eclipse/ 
^ available at http://www.na.astro.it/~radovich/wifix.htm 
available at http://www.ifa. hawaii.edu/~kaiser/imcat/l 
The usage of our pipeline can easily be learned within a 
few days by users having good experience in the reduction of 
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Table 1. Listed are the different UNIX architectures on 
which we used our pipeline so far. Besides the system we 
hst the compilers and the co-addition module used. We 
compared the final co-added images from different sys- 
tems with respect to flux and shape measurements from 
extracted objects. The individual reductions do not show 
significant differences. 



System 


compiler 


CO- addition 


Linux/i386 


gcc-2.95.3 


SWarp, EIS Drizzle 


Linux/ alpha 


gcc-2.96 


SWarp 


Linux/ AMD64 


gcc-3.3.4 


SWarp 


SunOS 


gcc-2.95.2 


SWarp, EIS Drizzle 


IBM AIX 


gcc-2.95.2 


SWarp 


HP-UX 


gcc-2.7.2.2 


SWarp 


OSFl 


gcc-2.7-97rl 


SWarp 


IRIX64 


gcc-3.0.4 


SWarp 



our tools complementary to the survey systems devel- 
oped by other groups such as TERAPIX^^ or ASTRO- 
WISE^'^. We note tha t the EI S team recently released a 
WFI pipeline l)Vandam&.2002i. l2004j) with a functionality 
similar to ours^'*. An alternative, and widely used reduc- 
tion package for CCD mosaics is the mscred environment 
within IRAF. It is described in iValdea (,2002.1 . 

2.3. Parallel processing 

A great deal of the reduction processes is performed in- 
dependently on individual chips or on a set of images be- 
longing to the same CCD. Only very few steps require 
information or imaging data from the complete mosaic. 
This mainly concerns the estimation of a global astromet- 
ric and photometric solution for the entire mosaic (which 
is done on catalogues) and the final image co-addition. 
Hence, most processing tasks can naturally be performed 
in parallel if the mosaic data is split up on a detector ba- 
sis and if a multi-processor architecture is at hand. Our 
pipeline is designed to perform this kind of parallelisation 
on a single multi-processor machine as well as on machines 
connected over a network (such as a Beowulf Cluster for 
instance). This parallelisation scheme of performing the 
same task simultaneously on different machine processors 
rather than using explicit parallel code for the individual 
tasks has two main advantages: 

— No new software has to be developed for the paral- 
lel processing. We developed a parallel manager that 
launches and surveys script execution on the individ- 
ual cluster nodes. The parallel manager is also used in 
standard, single-node processing so that the pipeline 

single-chip cameras, as was experienced by several guests who 
made use of our visitor programme. 

http://terapix.iap.fr/ 

http://www.astro-wise.org/ 

http://www.eso.org/eis/ 



Table 2. Processing times for a WFI@2.2m set consist- 
ing of 10 Bias, 12 Dark, 12 Skyflat and 11 Science frames. 
The third column lists the percentage of the processing 
time with respect to the single-node case. The number 
in parentheses gives the percentage if the whole process 
could be parallelised, i.e. 100% divided by the number 
of nodes. The fourth column lists the formal number of 
Science frame pixels processed in one second. The times in- 
clude the full processing from raw data as they arrive from 
the ESO archive to the final co-added images [it does not 
include the time for eye-balling and masking of images (see 
Sect. l3.in|l ]. We note that the parts that cannot be paral- 
lelised (an estimate for the global astrometric alignment 
with Astrometrix, the global photometric solution and 
the final co-addition with SWarp) and have to be done on 
a single processor add up to 24 min in our example. This 
becomes a very significant fraction of the total processing 
time in the 4 and 8 node cases. Each node is equipped 
with an Athlon XP 2800+ processor and 1 GB of main 
memory. For a comparison, for the current EIS pipeline a 
throughput of .4-0.5 Mpix/s o n a double-processor ma- 
chine is quoted l)Vandamell200l . 



used nodes 


time [min] 


time [%] 


Kpix/s 


1 


240 


100 (100) 


50 


2 


131 


54.0 (50) 


90 


4 


77 


32.0 (25) 


152 


8 


52 


21.7 (12.5) 


226 



operation is unified in single- and multi-processor 
mode. 

— As the operations on different chips are completely in- 
dependent and no communication is necessary between 
the tasks, no message passing/synchronisation or data 
exchange schemes have to be implemented. We only 
synchronise the processing after every processing step, 
waiting until a certain job is done on all chips before 
the next task is launched. 

The major disadvantage in our scheme not allowing a par- 
allelisation in the processing of a given chip is that re- 
sources are not used optimally if more nodes/processors 
than chips are available (nodes are simply staying idle for 
the complete process) or if the number of nodes is not 
an integer multiple of chips (some nodes have to work 
on more chips than others, basically doubling processing 
time). Nevertheless, we think that the advantages heavily 
outweigh this disadvantage for our purposes. Moreover, 
most multi-chip cameras so far have 2^ chips matching 
typical architectures of computer clusters. In Table [3 we 
compare processing times from single- and multi-node pro- 
cessing of an observing run. 

2.4. Terminology 

We use the following terminology throughout the paper: 

— CCD, chip - one of the detectors in a WFI 
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— exposure - a single WFI shot through the telescope 

— image ~ the part of the exposure that belongs to a 
particular CCD 

— frame - this term is used as a synonym for an image. 

— BIAS - an exposure/image with zero seconds exposure 
time 

— DARK - an exposure/image with non-zero exposure 
time, keeping the shutter closed 

— FLAT " an exposure/image of a uniformly illumi- 
nated area; this can be the telescope dome giv- 
ing DOMEFLATs or sky observations during evening 
and/or morning twilight providing SKYFLATs 

— SCIENCE - an exposure/image of the actual target, 
not a calibration image 

— SUPERFLAT - properly stacked SCIENCE data to 
extract large scale illuminations or fringe patterns. 

— STANDARD - an exposure/image of photometric 
standard stars 

— other terms written in CAPITAL letters denote addi- 
tional images or calibration frames. Their meaning will 
be clear within the context. 

— mosaic - this term is used as a synonym for exposure 
or for the final co-added image. 

— dithering - offsetting the telescope between the expo- 
sures 

— overlap - images from different CCDs and exposures 
overlap if the dithering between the exposures was 
large enough 

~ stack - a set of n images belonging to the same chip 
that have been combined pixel by pixel. 

— Names of software packages are written in the 
Typewriter font. 

3. Pre-reduction {Run processing) 

In the following we describe our algorithms for the pre- 
reduction of optical data, i.e. the removal of the instru- 
mental signature. The first issue is the compilation of data 
for this step. For most of the effects to be corrected for 
(instrument bias, bad CCD pixels, flat-field properties) 
we can safely assume stability of the instrument and the 
CCD characteristics over several days or even a few weeks. 
Hence, in many cases we can collect data from a complete 
observing run and benefit from having many images for 
the necessary corrections from which most are of statisti- 
cal nature. As described below, the matter can become 
more complex when dealing with strong fringes in red 
passbands. Here, the time scales from which SCIENCE 
data can safely be combined in the process is much shorter, 
sometimes only a couple of hours. The issue is further dis- 
cussed in Sect. 13.91 In any case we say that we perform 
the pre-reduction process on a -Run basis regardless of how 
long this period actually is. The pre-reduction is done in- 
dependently on each CCD of a mosaic camera. Only in one 
step, the sky-background equahsation (see Sect. 13.61) . the 
action to be performed on a CCD depends on properties 
from the rest of the mosaic. Hence, unless stated other- 
wise, each step described below has to be performed on a 



detector basis. For this part of the pipeline we can use two 
different software packages, one based on FLIPS, the other 
on Eclipse. FLIPS is very I/O intensive but has very low 
memory requirements, while Eclipse reduces the neces- 
sary I/O to a minimum and operates very efficiently on 
imaging data by keeping it in virtual memory. Depending 
on the size of the data set and the computer equipment at 
hand, one or the other is preferable. In the following we 
focus on the description of the Eclipse package and we 
mention FLIPS where its functionality differs significantly 
from Eclipse. 

3.1. Handling the data and the FITS headers 

The variety of ways in which header information and raw 
data from WFIs are stored in FITS files are as large as the 
number of instruments, i.e. so far there is no established 
standard on the FITS format of CCD mosaics. In order to 
cope with the different formats and to unify the treatment, 
we perform the following tasks on raw data: 

1. If the data are stored in Multiple Extension Fits 
(MEF) files, we extract all the images from them. All 
subsequent pipeline processing is done on individual 
images and also our pipeline parallelisation is based 
on the simultaneous execution of the same task on dif- 
ferent images. 

2. We substitute the original image headers from the 
chips by a new one containing only a minimum set of 
keywords necessary for a successful pipeline processing 
(see App. EJ- Especially the astrometric calibration 
depends on correct header entries. 

3. If necessary we flip and/or rotate individual CCDs to 
bring all images of the mosaic to the same orientation 
with respect to the sky. Only rotations by 90 degrees, 
that do not require pixel resampling, are performed. 

All these tasks are performed by a qf its-based utility. 

3.2. Modes, medians, and the stacking of images 

In most pre-reduction steps we have to perform robust 
statistical estimates on a set of pixels following Poisson 
statistics. At some points we have to estimate the sky- 
background and its variance in SCIENCE data. A typical 
distribution of pixel values on a SCIENCE image is shown 
in Fig. 121 The second common operation is to statistically 
stack several input images to produce a noise-reduced cal- 
ibration frame. This master image is then used for subse- 
quent calculations on SCIENCE images. In the following 
we define several operators acting on images: 

— value=median(image): 

The median of an image is estimated. This is done 
by collecting a representative sample from the CCDs 
pixels (we take about 5% of the detector pixels in a 
region around the CCD centre), sorting the obtained 
pixel array and returning its middle element (in case 
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Fig. 2. Shown is a typical pixel value distribution in an astro- 
nomical image. The clear peak denotes the brightness of the 
night sky. Due to the presence of objects the distribution is 
strongly skewed towards high values. The most accurate es- 
timate for the sky comes from analysing the histogram and 
determining the peak (the mode of the distribution) directly. 
The median (for a set of A'' points the median element is de- 
fined so that there are as many greater as smaller values than 
the median in the set) also gives a very robust (and for our 
purposes sufficient though biased) estimate for it. Due to the 
long tail of the distribution, a straight mean is completely use- 
less as an estimate for the sky value. As the values to the left 
of the mode form half of a Gaussian distribution, an estimate 
for the sky variance can be obtained from that part. 

the array has an even number of elements we return 
the mean of the two middle elements). 

— value— mode(image): 

The mode of an image is estimated. We consider the 
same representative pixel sub-sample as for the median 
estimation, build a smoothed histogram and return the 
peak value. 

— image=rescale (images ) : 

This operator rescales a set of images so that they have 
the same median after the process. The resulting me- 
dian is chosen to be the mean of the medians (we write 
meanmed for it) from the input images. Hence on each 
input image we perform the operation image image 
* meanmed / median(image). The operator is usually 
applied to a set of images before they are stacked (as 
SKYFLATs or SCIENCE images). The median of the 
stacked image is then by construction equal to mean- 
med. 

— image— stack(images): 

A set of input images is stacked to produce a noise- 
reduced master image. The following procedure is per- 
formed independently on each pixel position. We col- 



lect all pixel values from the input images into an array, 
sort it and reject several low and high values (typically 
we reject the three lowest and three highest values if 
we have 15 input images or more). From the rest we 
estimate the median that goes into the master image. 
Here FLIPS uses a more sophisticated algorithm. On 
the remaining array it first performs an iterative sigma 
clipping. It estimates mean and sigma, rejects low and 
high values (typically pixels lying more than 3.5 sigma 
below and above the mean) and repeats the procedure 
until no more pixels are rejected. From the rest it re- 
turns the median for the master image. 

We will use additional pseudo operators in the following 
whose meaning and behaviour will be clear within the con- 
text. All calculations with these operators are written in 
slanted notation. 

3.3. A first quality check, Overscan correction, master 
BIASes and DARKs 

Before any exposure enters the reduction process, we es- 
timate its mode. If this estimate lies outside predefined 
boundaries, the exposure will not be processed. This re- 
jects most of the completely faulty images (such as satu- 
rated SKYFLATs) at the very beginning. After this ini- 
tial quality check, the first step in the reduction process 
of each image is the correction of an overall estimate for 
the BIAS value by considering pixels in not illuminated 
parts of each CCD (the overscan region). This first-order 
BIAS correction is done by collecting for each line all pix- 
els in the overscan region, rejecting the lowest and high- 
est values (usually the overscan regions contains about 40 
columns and we reject the 5 lowest and 5 highest values) 
and estimating a straight mean from the rest. This mean is 
subtracted from the corresponding line. After this correc- 
tion, the overscan regions are trimmed off the images. For 
correcting spatially varying BIAS patterns in the FLAT 
or SCIENCE images, a master BIAS is created for each 
CCD from several individual BIAS exposures: 

1. Each BIAS exposure is overscan corrected and 
trimmed. 

2. The master BIAS is formed by masBIAS^stack(BIAS 
frames ) 

In the same way, a master DARK (masDARK) frame 
is created. So far, we do not correct FLAT or SCIENCE 
frames for a possible dark current but we use the mas- 
DARK for the identification of bad pixels, rows and 
columns (see Sect. 13)). 

3.4. Methods for flat-fielding 

The illumination of a CCD and its sensitivity is in gen- 
eral not homogeneous over the detector area. Illumination 
usually varies on large scales, but sensitivity can change 
significantly from pixel to pixel (see Fig. Ej). In order to 
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Fig. 3. A skyflat from the CFH12K@CFHT camera dCuillandre et aLlf20n(fl . We note the significant intrinsic sensitivity and 
gain differences between the chips. The counts of the second left detector in the top row are about a factor 1.3 higher than that 
of the two chips at the corners in the bottom row. Within a CCD, large-scale gradients are up to 8%. Especially the four chips 
on the right also show small-scale sensitivity variations across the whole CCD. 



overcome these effects and to achieve a homogeneous sur- 
face brightness, SCIENCE images have to be flat-fielded. 
The flat-field pattern typically varies with the observa- 
tion configuration and hence calibration images have to 
be obtained for each filter to be used. There are three 
major possibilities to obtain an estimate for the flat-field- 
pattcrn: 

1. DOMEFLATs: 

One can obtain flat images by pointing the telescope 
to an illuminated, blank surface in the dome. These 
DOMEFLATs can be taken during day time, the im- 
ages do not contain any astronomical sources, we can 
take as many exposures as we wish and we can ad- 
just the count rate in the frames. This ensures FLATs 
with a very high signal-to-noise. However, especially 
for WFIs, the uniform illumination of a large enough 
area is very difficult to achieve. In addition, lamps emit 
light at a maximum temperature of around 3000 K so 
that the spectral distribution is very red. In red fil- 
ters this can lead to an excess amount of fringing as 
compared to SKYFLATs. For these reasons, the flat 
calibration of SCIENCE frames with DOMEFLATs is 
usually inferior to the two methods described below 



(especially at the correction of large-scale illumination 
variations). We only use them if SKYFLATs are not 
available. 
SKYFLATs: 

One can also obtain flat images by observing the sky 
with short exposure times during the evening and 
morning twilight. These flats resemble the night sky 
during SCIENCE observations better than dome flats, 
although for WFIs covering an area of 30' x 30' or 
more we expect variations of the twilight sky over their 
field-of-view. At many telescopes, an automatic setup 
to obtain exposures with a desired count-rate are avail- 
able and hence thes e fiats also have high S/N. See also 
iTvson fc C^l^^993^ for a guide to obtain twilight flats 
if an automatic procedure is not available. However, 
only a very limited number can be obtained during 
the twilight phase. This holds especially if several fil- 
ters have to be calibrated and/or if the readout time of 
the detector is rather long. Using an exposure time too 
short in bright twilight can lead to a non-uniform illu- 
mination due to shutter effects; see Sect. 13.91 However, 
CCDs are usually very stable over a few days or weeks 
so that a larger number of SKYFLATs taken during 
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several nights can be combined into good master flat 
images. 
3. SUPERFLATs: 

In addition, one can try to extract the flat-field pat- 
tern from the SCIENCE observations itself. If a suffi- 
cient number of SCIENCE exposures is at hand (usu- 
ally more than a dozen), and if the dither pattern was 
significantly larger than the largest object in the field, 
each pixel of the camera will see the sky-background 
several times. Hence, a proper combination of these ex- 
posures yields a master FLAT that closely represents 
the night sky during observations. We refer to such a 
FLAT as a SUPERFLAT. A straightforward applica- 
tion of this method is often hampered by several fac- 
tors: During phases of grey and bright time, the night 
sky shows large gradients and variable reflections in the 
dome and telescope can occur. Thus, a careful selec- 
tion of images that go into the SUPERFLAT has to be 
done. In medium/narrow band filters and in the ultra 
violet, the counts and hence the S/N of SUPERFLATs 
in these bands are typically low. Furthermore, the tech- 
nique is very difficult to apply in programmes where 
the imaged target has a size comparable to the field- 
of-view. In this case also large dither patterns cannot 
assure sufficient sky coverage on the complete mosaic. 

As the superflat technique cannot be applied in general 
we use the following two-stage flat-fielding process: 

1. In a flrst step, the SCIENCE observations are flat- 
fielded with SKYFLATs (if those are not at hand we 
use DOMEFLATs instead). This typically corrects the 
small-scale sensitivity variations and leaves large-scale 
gradients around the 3% level on the scale of a chip 
(see Fig. El . 

2. If the data permit it, a SUPERFLAT is created out of 
the flat-fielded SCIENCE observations, smoothed with 
a large kernel to extract remaining large-scale varia- 
tions and is applied to the individual images. For our 
empty field observations at the WFI@2.2m this leaves 
typical large-scale variations around the 1% level in 
U,B,V,R broad band observations (in / the presence 
of strong fringing often leads to significantly higher 
residuals). 

This two-stage flat- fielding process is v ery similar to the 
method adopted bv lAlcala et al 

3.5. The creation of master DOME-/SKYFLATs 

The master FLAT (SKYFLAT or DOMEFLAT) for each 
CCD is created as follows: 

1. All individual FLAT exposures are overscan-corrected 
and trimmed. 

2. The masBIAS is subtracted from all images. 

3. The FLATs are rescaled to a common median: 
rescFLAT=rescale(FLAT). 

4. We form the master FLAT by mas- 
FLAT=stacJf('rescFLAT frames) 



3.6. Sky-background equalisation 

Within a mosaic, CCDs have different quantum efficien- 
cies, varying intrinsic gains and hence different photomet- 
ric zero points (see Fig. . For the later photometric cal- 
ibration it is desirable to equalise the photometric zero 
points in all detectors, which is achieved by scaling all 
CCDs to the same sky-background. We rescale all chips to 
the median of the CCD with the highest count-rate dur- 
ing flat-fielding. If possible we perform this step within 
the SUPERFLAT correction as its median estimation is 
more robust than with the SKYFLATs. This is because 
SKYFLATs show larger variations in brightness than the 
SUPERFLATs which are calculated from already flat- 
fielded data. We estimate that photometric zero points 
of the mosaic agree with an rms scatter of about 0.01-0.03 
mag after sky-background equahsation (see Sect. l4.6|l . 

3.7. SCIENCE image processing, the creation and 
application of the SUPERFLAT 

The individual steps in the SCIENCE frame processing 
are: 

1. The images are overscan-corrected and trimmed. 
Afterwards the masBIAS is subtracted and the 
frames are divided by masFLAT. We write 
aatSCIENCE= (SCIENCE- masBIAS ) /masFLAT. 

In the case where no SUPERFLAT correction is 
applied to the data, masFLAT is rescaled at this step 
to equalise the sky-background between individual 
detectors (see Sect. I3.6|l . 

2. For the creation of the SUPERFLAT we first remove 
astronomical sources from the flat SCIENCE images. 
To th is end we run SExtractor ijBertin fc Arnoutl 
Il996|) on them and create a new set of im- 
ages where pixels belonging to objects (i.e. 
above certain detection thresholds) are flagged 
{objSCIENCE=BatSCIENCE- OBJECTS). See also 
Fig. 0] The flagged pixels are not taken into account 
in the subsequent processing. 

3. The obj SCIENCE images are rescaled to a common 
median rescobjSCIENCE=rescale(objSCIENCE). 

4. We calculate the SUPERFLAT 
[SUPERFLAT^stack(rescobjSCIENCE frames )] . 
If all input pixels of a given position are flagged (i.e. 
all the images had an object at that position) we 
assign the meanmed value from the obj SCIENCE 
frames to the SUPERFLAT. 

5. The SUPERFLAT is heavily smoothed by creating 
a SExtractor BACKGROUND che ck-image with a 
background mesh of 512 pixels fsee iBertinI l2003l on 
the SExtractor sky-background estimation). This im- 
age, the illumination correction image, forms the 
basis for removing large-scale flat-field variations 
(ILL UMINATION^smooth (SUPERFLAT)) . 

6. The flatSCIENCE images are divided by the 
ILLUMINATION image which has been rescaled 
to equalise the sky-background of the differ- 
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ent detectors (see Sect. We write illum- 

SCIENCE=RatSCIENCE/rescaled(ILLUMINATION). 

For blue passbands not showing any fringing, the pre- 
reduction ends here and the illumSCIENCE images are 
used in the subsequent processing. See Fig. [S] for an ex- 
ample of a pre-reduced F-band exposure. 

3.8. Fringing correction in red passbands 

Fringing is observed as an additional, additive instrumen- 
tal signature in red passbands. It is most prominent on 
cameras that use thinned CCDs and hence are optimised 
for observations in blue passbands. Fringes show up as spa- 
tially quickly varying, wave-like structures on the CCDs 
(see Fig. EJ. The geometry of these patterns usually does 
not change with time since this interference effect is cre- 
ated in the substrate of the CCD itself. WFI@2.2m shows 
fringes with an amplitude of about 1% as compared to 
the sky-background in broad-band R. In the /- and Z- 
bands, fringing becomes much stronger and reaches up to 
about 10% of the night sky. Unfortunately, contrary to 
the geometry of the fringing, its amplitude is not stable 
in consecutive SCIENCE exposures since it strongly de- 
pends on the night sky conditions (as brightness, cloud 
coverage), the position on the sky and the airmass. If a 
reasonably good SUPERFLAT can be constructed from 
a sufficient number of SCIENCE frames obtained under 
stable sky conditions, a possible way to correct for fringes 
is the following: 

1. Besides the large-scale sky variations not corrected for 
by SKYFLATs, the SUPERFLAT contains the fringe 
pattern as an additive component. Hence, the fringe 
pattern can be isolated from the SUPERFLAT by 
FRINGE^SUPERFLAT~ILL UMINATION. 

2. Individual SCIENCE frames are corrected for 
fringing by fringeSCIENCE=illumSCIENCE-f x 
FRINGE. We assume that the fringe amplitude di- 
rectly scales with sky brightness and f is calculated by 
f=median (illumSCIENCE ) /median (ILL UMINATION) ^ 

In the case of WFI(Q)2.2m, this method removes very ef- 
ficiently the low-level fringes in the i?-band. Fringing is 
usually no longer visible and we estimate possible residu- 
als well below 0.1%. In the case of /- and Z-bands, fringes 
can be suppressed to a level of about 0.1% of the sky am- 
plitude if a very good SUPERFLAT can be constructed 
(see Fig.|S}. If this is not the case, our approach may per- 
form very poorly in reducing the fringe amplitude and we 
cannot propose a pipeline solution to the problem at this 
stage. 

3.9. Guidelines for constructing calibration images 

The success of the image pre-reduction heavily depends 
on the quality of calibration data at hand. In the creation 
of flatSCIENCE, we propagate the noise in masBIAS and 
masFLAT to our SCIENCE frames. As the noise in these 



calibration frames is of statistical nature we can dimin- 
ish it by using as many images in the stacking process as 
possible. For the successful creation of a SUPERFLAT, 
and hence for the later quality of large-scale illumination 
and fringing correction, not only the number of images is 
important but it is essential that each pixel in the mosaic 
sees blank sky several times during SCIENCE observa- 
tions. This suggests that the best observing strategy to 
achieve this goal is a large dither pattern between con- 
secutive exposures (ideally it should be wider than the 
largest object in the field). In the case that the target 
occupies a significant fraction of the mosaic (as big galax- 
ies or globular clusters for instance) the best strategy is 
to observe a neighbouring blank field for a short while if 
a SUPERFLAT and/or fringe correction is important. In 
the following we give some additional guidelines leading 
to good results in most of our own reductions: 

— For the stacking process for masBIAS (masDARK) and 
masFLAT, about 15-20 individual images of each type 
should be acquired. As the noise in the final images is 
inversely proportional to the square root of the number 
of input images this is a good compromise between 
the number of exposures that can be taken during a 
typical observing run (in the case of SKYFLATs) and 
the noise reduction. 

— The observer should find out the minimum exposure 
time required for the FLATs in order to suppress shut- 
ter effects. This can be done by taking DOMEFLATs 
or SKYFLATs with varying exposure time. One then 
normalises the flat fields and divides them. Any sys- 
tematic residuals are then due to shutter effects, and 
one can determine the minimum usable exposure time. 
This is illustrated in Fig. |7|for the prime focus imager 
of the 4.2m WiUiam Herschel Telescope (WHT) at La 
Palma, which has a comparatively slow iris like shutter. 
Thus, images with short exposure times receive more 
light in the centre than in the corners. In any case, 
a SUPERFLAT of the data would recover any such 
shutter-related inhomogeneities in the FLAT, but we 
recommend to get everything correct from the begin- 
ning if possible. 

For newer or future WFIs such as OmegaCAM, shut- 
ters are constructed in which two blades move with 
constant speed from one side to the other exposing 
and covering each pixel of t he detector arra y for ex- 
actly the same time (see e.g. lR,eif et al.ll200^ . 

— As mentioned above, the creation of SUPERFLATs 
is not trivial in general. This is primarily an issue of 
which SCIENCE images to include. In the case of blue 
passbands, where only the smoothed ILLUMINATION 
image is used to correct for large-scale illumination 
gradients, it is safe to construct it from all SCIENCE 
frames of several, consecutive nights. Blue passbands 
are typically observed during periods of dark time, 
and hence the sky conditions are sufficiently stable. 
Also, with the moon down we observed no significant 
dependence of large-scale sky variations on telescope 
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Fig. 4. For the creation of SUPERFLATs we mask pixels belonging to astronomical sources in the SCIENCE frames. The left 
panel shows a SCIENCE image before and the right panel after objects have been detected by SExtractor and flagged. We 
mask all structures having 50 contiguous pixels with 1.0 sigma above the sky-background. This helps to pick up also pixels of 
extended halos around bright stars. For images containing strong fringing we change the parameters to 7 pixels above 5 sigma 
as otherwise fringes would be detected and masked. 




Fig. 7. Ratio of two normalised sky flats taken with the 
prime focus imager of the WHT (2 CCDs). The images 
were exposed for 0.3 and 1.0 seconds, respectively. One 
can see the way the iris shutter exposed and covered the 
detector again. The illumination difference between the 
brightest and faintest part in this representation is 20%. 
For the WHT prime focus imager it is recommended to 
expose for at least 2.0 seconds. 



pointing. Hence, a robust SUPERFLAT can be ob- 
tained with a sufficient number of images during an 
observing run. This is no longer the case if we need 
to correct for strong fringing. We observe that our as- 



sumption of a sole dependence of the fringe amplitude 
on sky brightness typically fails when constructing the 
SUPERFLAT from SCIENCE frames of different sky 
positions. As the appearance of the overall pattern is 
very stable, this means either that the fringe ampli- 
tude is no longer a linear function of sky brightness or 
that the sky-background varies on scales significantly 
smaller than a CCD so that the scaling factor becomes 
dependent on CCD position. In addition, a similar be- 
haviour is observed as a function of time, depending 
on the excitation of the 0H~ night-sky emission lines. 
For our blank-field observations we obtain good results 
in the fringing correction if we observe the same target 
between 10 and 15 times with a large dither pattern 
within an hour (at WFI@2.2m with an overhead of 
about 2 minutes per image this can be achieved with 
300 s exposures). 

If a very good SUPERFLAT could be constructed, flat- 
fielding results from our proposed method and from a 
pure SUPERFLAT application (i. e. using flat-fielding 
method|3|in Sect. 13. 4|) should be compared if flatness in 
SCIENCE images is crucial. Especially if only a small 
number of individual frames went into the construction 
of masFLAT, the direct SUPERFLAT approach often 
gives better results. 

Offsetting the telescope between the exposures is fun- 
damental for high-quality, high-S/N mosaics. The 
dither box, i.e. the box that encloses all dither off- 
sets, should be clearly larger than the gaps between 
the CCDs and the objects in the field. This has sev- 
eral advantages as compared to the staring mode (no 
offsets at all) or to the application of only small offsets. 
First, the CCDs in a multi-chip camera are fully in- 
dependent from each other. They see the sky through 
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Fig. 5. The pre-reduction steps on a l/-band exposure from WFI@2.2m: Panel (A) shows the raw image, panel (B) the result 
from applying masBIAS, trimming the CCDs and flat-fielding with the masFLAT. The high S/N masFLAT takes out small-scale 
variations but leaves large-scale residuals of up to 3% on the scale of the CCD [panel (D)]. These and the differences in the sky 
count-rate are removed after the application of ILLUMINATION giving a flatness of 1% over the entire mosaic in most of the 
cases as seen in panels (C) and (E). 



difFerent sections of the filter, and they have their in- 
dividual flat-fields, gains and read-out noise. Choosing 
a wide dither pattern is the easiest way to establish an 
accurate global photometric and astrometric solution 
for the entire mosaic, based on enough overlap objects. 
Second, the wide dither pattern allows for a signif- 
icantly better superflattening of the data, since the 
objects do not fall on top of themselves in the stacks 



and thus for every pixel a good estimate of the back- 
ground can be obtained. Besides, remaining very low- 
amplitude patterns in the sky background caused by 
improper flat-fields etc. do not add in the mosaic, but 
are averaged out. Thus, a wide dither pattern will lead 
to an improved sky background from which the S/N 
will benefit; see also Fig. [S] 
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Fig. 6. Fringing correction in WFI@2.2m 7-band data: For the shown correction, fifteen 300 s exposures have been obtained 
within 1.5 hours during stable, photometric conditions with a large dither pattern (3 arcmin in Ra and Dec) on a blank-field 
target. The upper left panel shows a SCIENCE frame before fringing correction, the upper right panel the extracted fringe 
pattern from the SUPERFLAT and the lower left panel the SCIENCE image after the correction. The lower right panel shows 
a zoom of a representative region (marked in the left two panels) after the correction. The residual fringes have an amplitude 
of about 0.1% of the night sky. 



Third, the object S/N in a mosaic with a (wide) dither 
pattern is significantly superior to one obtained in 
staring mode. The reason for this is that the mas- 
ter BIAS and the master FLAT are not noise-free, 
since they are created from a finite number of im- 
ages. An identical copy of this noise is then created 
in each SCIENCE image during the pre-processing. If 
no dither offsets are apphed, then the calibration noise 
in the N SCIENCE images is stacked on top of itself 



during the co-addition, and thus increases in the same 
way as the flux of the objects (oc N) instead of aver- 
aging out (oc ^/N). This is enforced by the use of an 
auto-guider, which keeps the sources in subsequent im- 
ages exactly at the same pixel position. Depending on 
the ratio of the calibration noise and the noise in the 
uncalibrated images themselves, the effective exposure 
time of the mosaic can be very significantly reduced. 
This is especially true for observations gained under 
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Fig. 8. Sky-background variations in a mosaic obtained from observations of several fields in staring mode (left), 
and for one obtained with a wide dither pattern (right). The contrast scaling for both mosaics is the same. Thus, 
science applications that require a very uniform background, e.g. the search for low surface brightness galaxies, profit 
in general substantially from large dither patterns. 



excellent (dark) conditions, since the sky noise in the 
SCIENCE images is then reduced, and the noise from 
the calibration frames becomes more dominant. One 
should therefore always aim at a sufficient number of 
calibration exposures, and go for larger dither offsets. 
This holds for single-chip cameras, too. 



3.10. Image Eyeballing and Manual Masking 

After the pre-reduction we visually inspect all SCIENCE 
images to identify potential problems needing manual in- 
tervention: 



— Faulty exposures (e.g. data taken out of focus) having 
passed the initial quality check fSect. rOjl are removed 
at this stage. 

— Problems of the pre-reduction process are identi- 
fied and necessary steps are iterated (this concerns 
mostly the selection of SCIENCE frames entering the 
SUPERFLAT). 

— The most important aspect of this manual pass 
through all SCIENCE data is the identification of 
large-scale image defects that must not enter the fi- 
nal co-addition process and which have to be manually 
masked at this stage. This will be described in more 
detail in Sect.El See also Fig. El 



3.11. Standard star processing and absolute 
photometric calibration 

Besides the SCIENCE images we have to process 
STANDARD frames to obtain an absolute photometric 
calibration. A fully automated absolute photometric cal- 
ibration up to the final co-added image is not yet imple- 
mented in our system. As a first step in this direction, 
the pipeline allows the photometric calibration of indi- 
vidual nights of a Run. In this section we describe how 
STANDARD exposures are pre-reduced and analysed to 
extract photometric parameters. If not stated otherwise, 
all magnitudes in this paper are given in the Vega system. 

STANDARD exposures are reduced in a fashion sim- 
ilar to SCIENCE exposures. They are first overscan cor- 
rected, then the masBIAS is subtracted. Flat-fielding is 
done with the masFLAT generated in the pre-reduction of 
the SCIENCE exposures and the SUPERFLAT image cre- 
ated from the SCIENCE exposures is applied to the pre- 
reduced STANDARD exposures. Using the SUPERFLATs 
created from the SCIENCE images is important to en- 
sure photometric homogeneity between SCIENCE and 
STANDARD exposures. 

If necessary the FRINGE maps created from the 
SCIENCE images are rescaled and subtracted from the 
superflatted STANDARD exposures. This step is particu- 
larly difficult and often leads to non-satisfactory results in 
the /-band. The reason is that STANDARD exposures are 
usually taken at different times during the night (some- 
times even in twilight), at different sky positions and 
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Fig. 9. Examples of large-scale defects: The upper right panel shows a satellite track, the other three panels various 
forms of reflections from bright stars inside and outside a WFI@2.2m field. All such features are typically unique to 
individual frames and do not appear in subsequent, sufficiently dithered exposures. Currently we have to register and 
to mask these defects by hand to prevent that they enter the final image co-addition (see Sect.O. Masks are built 
using the polygon option of the d.s9 visualisation tool. 



various airmasses. This often leads to fringing patterns 
with intensities very different from those present in the 
SCIENCE exposures. 

To generate object catalogues that can be used for pho- 
tometric calibration bad pixel masks for the STANDARD 
images are created. They mark pixels not suitable for pho- 
tometric measurements, such as bad pixels or cosmic rays. 
The creation of these masks (FLAG images henceforth) is 
described in more detail in Sect. El 

Using the FLAGs object catalogues are created with 
SExtractor from the reduced STANDARD ex pos ures. To 
match the found objects with iLandoltl l)l99 2') or 'St etsonI 
l)2000|) standard star catalogues an astrometric solution 
for the STANDARD images has to be derived first. This is 
done with the LDAC astrom tools to match objects found 
in the STANDARD images to the USN0-A2 catalogue. 



The number of objects found in the STANDARD expo- 
sures is sufficiently high to derive a second-order astro- 
metric solution separately for each image. This is enough 
to reliably merge standard star catalogues with our cata- 
logues. All fluxes in our STANDARD fields are normalised 
to an exposure time of 1 s. 

From the database of matched standard stars, the co- 
efficients of the equation, 

Mag - Magin.t = ZP + CE * airmass + CT * CI , (1) 

i.e. the zero point ZP, the extinction coefficient CE, and 
the colour term CT are determined, where Mag is the 
standard star's magnitude, Magj^g^. is the instrumental 
magnitude measured on the reduced standard frame (we 
use the SExtractor MAG_AUTO parameter for this pur- 
pose), and CI is a colour index, e.g. {B — V). This is done 
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using a non-linear least-squares Marquardt-Levenberg al- 
gorithm with an iterative 3a rejection, which allows re- 
jected points to be used again if they are compatible with 
later solutions. As this algorithm is not guaranteed to con- 
verge, the iteration is aborted as soon as one of the fol- 
lowing three criteria is true: 

1. The iteration converged and no new points are re- 
jected. 

2. The maximum number of iterations (set to 20) is 
reached. 

3. More than a fixed percentage (set to 50%) of all points 
are rejected. 

In a first step all three coefficients are fit simultaneously. 
However, in order to reliably estimate the extinction co- 
efficient standard star observations must be spread over 
a range of airmasses. This is sometimes neglected by ob- 
servers. To find an acceptable photometric solution in this 
case, the user can supply a default value for the extinction 
coefficient. The fit is then repeated with the extinction co- 
efficient fixed at the user supplied value and the zero point 
and colour term as free parameters. 

Although wide-field observations of Landolt/Stctson 
fields typically cover a wide range of stellar colours, the 
user can also supply a default value for the colour term. 
This is then used for a 1-parameter fit in which the zero 
point is the only free parameter. 

In an interactive step the user can then choose between 
the 1-, 2-, and 3-parameter solution, or reject the night as 
non-photometric. The FITS headers of the frames belong- 
ing to the same night are then updated with the selected 
zero point and extinction coefficient or left at the default 
value —1.0, indicating that no photometric calibration for 
that frame is available. 

We note that we perform the fit simultaneously for 
all mosaic chips. As discussed in Sect. 14.61 zero points 
of individual chips agree within 0.01-0.03 mag after sky- 
background equalisation. We do not take into account 
possible colour term variations between individual images 
that are expected due to slightly different CCD transmis- 
sion curves. Notable differences can occur in the U and Z 
filters that are cut off by the CCD transmission. 

Fig. llOlillustrates our photometric calibration which is 
implemented in Python. We compared our own zero point 
and extinction measurements for several nights with those 
recently obtained by the EIS team. TablclHlshows that the 
measurements are in very good agreement. 

3.12. Runs and Sets 

As discussed in the beginning of Sect. |31 the pre-reduction 
is done on a Run basis usually containing observations 
from different patches of the sky. Before proceeding we 
need to split up the SCIENCE exposures of a Run into the 
Sets that have to be co-added later. By Set we mean the 
series of all exposures of the same target in a particular 
filter. This means that all the following reduction steps 
up to the final co-addition have to be done on each Set 



independently. We note that a Set can have been observed 
in multiple Runs so that all Runs containing images of a 
certain Set have to be processed at this stage. 

4. Astrometric and photometric Set calibration 

4.1. Astrometry 

After the pre-processing, a global astrometric solution and 
a global relative photometric solution is calculated for all 
SCIENCE images. This is where the reduction of WFI 
data becomes much more complicated than the one for 
single-chip cameras. 

In a first step, high S/N objects in each image are de- 
tected by SExtractor, and a catalogue of non-saturated 
stars is generated. Typically, we use all objects hav- 
ing at least 5 contiguous pixels with 5a above the sky- 
background in the following analysis (these numbers may 
vary according to filter and exposure time; in the ?7-band 
for instance, we need to lower these thresholds in order to 
have enough sources to compare with a standard star cat- 
alogue). This usually gives us between 3 and 6 objects per 
square arcmin in high-galactic latitude empty field obser- 
vations. Based on a comparison with the USN0-A2 (see 
iMonet et al .119981) astrometric reference catalogue (or any 
other reference catalogue), a zero-order, single shift astro- 
metric solution is calculated for each image. For a single- 
chip camera with a small field-of-view such an approach is 
often sufficient, but it no longer holds for multi-chip cam- 
eras with their large field-of-view. In this case, the CCDs 
can be rotated with respect to each other, tilted against 
the focal plane, and in general cover areas at a distance 
from the optical axis where field distortions become large. 
Fig. shows the difference between a zero order (sin- 
gle shift with respect to a reference catalogue) and a full 
astrometric second-order astrometric solution per image. 
From this figure it is obvious that the simple shift-and-add 
approach will not work for the entire mosaic. The issue is 
further complicated by the gaps between the CCDs and 
large dither patterns that are used to cover them. Thus, 
images with very different distortions overlap. In addition, 
due to the large field-of-view, one must take into account 
that the spherically curved sky is mapped into the fiat 
tangential detector plane. 

In the second step, Mario Radovich's Astrometrix^^ 
package is used to determine third-order polynomials for 
the astrometric solution of each image. This corrects for 
the aforementioned effects, and thus allows for the proper 
mosaicing in the later co-addition process. For this pur- 
pose all high S/N objects (stars and galaxies) detected in 
the first step are matched, including those from the over- 
lap regions. The latter ones are most important in estab- 
lishing a global astrometric (and photometric) solution, 
since the accuracy of available reference catalogues such 
as the USN0-A2 with an rms of about O'.'S is insufficient 
for sub-pixel registration. Thus the astrometric solution 

http://www.na.astro.it/~radovich/WIFIX/astrom.ps 
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Fig. 10. The plots illustrate our zero point estimation in broad-band B during a clear night. Eleven standard field observations 
with a good airmass coverage were obtained giving a total of about 4000 individual standard star measurements (the Stetson 
standard star catalogue was used here). Points indicate the measurements, the lines the result of a three parameter fit of eq. 

Table 3. The table shows comparisons of our own photometric zero point and extinction measurements with 
those from the EIS team. We independently calibrated several nights in the WFI@2.2m broad-band B fil- 
ter (ESO Filter ID: BB#B/123_ES0878). The EIS values were taken from the EIS XMM product logs (see 
[htt p:/ /www. eso.org/science/eis/surveys/releaseJ35000024_XMM. html) of the fields listed in column 2. From the 
comparison we excluded the measurement of one non-photometric night (30/06/2003). The 7th column lists 
AZP=Bonn_ZP-EIS_ZP and the 8th column the number of parameters fitted in eq. Values marked with an 
asterisk were not fitted but default values have been used instead. We note that the results are in very good agree- 
ment. 



Date 


Field 


EIS.ZP 


EIS ext. 


Bonn_ZP 


Bonn ext. 


AZP 


#par (Bonn) 


22/06/2003 


BPM16274 


24.80 


-0.199 


24.80 


-0.16 


0.00 


3 


05/03/2003 


HE1104-1805 


24.86 


-0.242 


24.91 


-0.30 


0.05 


3 


30/03/2003 


NGC4666 


24.71 


-0.184 


24.75 


-0.14 


0.04 


3 


06/08/2003 


LBQS^212-1759 


24.82 


-0.205 


24.88 


-0.22(*) 


0.06 


2 


01/04/2003 


Q1246-057 


24.64 


-0.122 


24.67 


-0.08 


0.03 


3 


27/01/2003 


SHARC-2 


24.58 


-0.108 


24.51 


-0.04 


-0.07 


3 



is determined from the data itself instead of the reference 
catalogue, which is usually based on data taken with a 
much smaller angular resolution than the images that are 
processed. The reference catalogue is used only to fix the 
solution with respect to absolute sky coordinates within 
its stated rms. A wide dither pattern is required to com- 
pute a reliable global astrometric solution, so that differ- 
ent regions of the CCD are mapped by overlapping sources 
and the astrometric solution is properly constrained. For 
a more detailed desc ription of the Astromet rix tool see 
iRadovichI ll2n02^ and lMcCracken et all t()n?t\ . 



4.2. Accuracy of our Astrometric calibration 

The astrometric calibration ensures that we can produce a 
deep mosaic out of individual images. The accuracy with 
which the images have to be aligned strongly depends on 
the scientific application in mind. In most cases, the final 
goal is the measurement of moments M™" from the light 
distribution I{6) of objects, that is: 

7^"":= J d^e9'P91^W{e)I{e);m,n€Afo 

where W{9) is some weighting function, 9 = \6\ and 
{9i,92) ~ (0,0) is the centre of the object, which is cho- 
sen to be the point where the weighted dipole moments 
are zero. In our weak lensing studies, we have to estimate 
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Fig. 11. Difference in object position between a single-shift 
approach and a full two-dimensional second-order astrometric 
solution for the WFI@2.2m. In other words, shown are the 
higher-order terms needed for matching the images to the sky. 
The patterns belonging to the left two CCDs are due to a 
rotation with respect to the mosaic. The maximum position 
difference in the plot is about six pixels, still a fairly small 
value compared to other telescope designs. It becomes clear 
that a single, global distortion polynomial for all 8 CCDs does 
not work. Instead, every image has to be treated individually. 



accurately moments up to the fourth order. The shape of 
an object is constructed out of the second moments qij, 
i.e.; 



where we choose W{9) = l/(27rr2) exp(-|6'| V(2r2)) with 
fg = 3 image pixels henceforth. The following two fourth- 
order moments are important for correcting object shapes 
for P SF effects in weak lensing studies (see iKaiser et al.l 
Il995l) : 



2We^ + 2W'r]l 2W'r]iTj2 
2W'r]iri2 2W9^ + 2W'r]j 



m 



and 



where r/i := {ej - 91), 772 := 29i92, 773 + 2W'9'^ 

and W' denotes the derivative with respect to 9^ . As the 
quantities under consideration are linear in the hght dis- 
tribution I{6), the measurement in the mosaic can be pre- 
dicted from the measurements in the individual images (if 
a linear image co-addition is performed). If the mosaic is 
constructed by a straight mean of N input images, the 



value 
equal to: 



measured in the co-added image has to be 



1 ^ 



fc=l 



where the q^j- are the moments in the individual images. A 
systematic deviation of the quantity qf-^ := qfj°^^^'^ — qf-^'^'^ 
(and in corresponding expressions for X^^ and X^'°^) from 
zero is hence a good indicator that the astrometric solu- 
tion should be improved. Results of this test for a deep 
i?-band mosaic constructed out of 30 exposures with sub- 
arcsecond seeing are shown in Figs. El ^1 and E| The 
agreement between the predictions and the measurements 
is in general very good, only the results of X^^ may indi- 
cate a small systematic offset in the co-added image. We 
conclude that the alignment accuracy of Astrometrix is 
sufficient for our weak lensing studies. For our WFI@2.2m 
i?-band observations Astrometrix formally estimates an 
internal astrometric accuracy of 1/10— 1/5 of a WFI pixel 
(0'.'02 - 0'.'04) over the whole field-of-view. This result 
is con sistent with the accuracy obtained bv lAlcala et alJ 
l|2002() for WFI@2.2m observations. However, this accu- 
racy strongly depends on filter, instrument and the ob- 
serving strategy, i.e. the number of bright stars that can be 
extracted, th e kind of distort i ons an d the dither pattern. 
For instance, iRadovich et alJ l|2004^ obtain an alignment 
accuracy of abou t 1/3 of a pixel for WFI @2.2m ?7-band 
observations and iMcCracken et al.l l)2003|) 2/5 of a pixel 
for B, V, R, I observations with the CFH12K@CFHT in- 
strument. Finally we compare the size of the PSF in the 
individual images and the mosaic. Fig. 1 151 shows that also 
the size of the PSF in the mosaic can very well be predicted 
by measurements in the individual exposures entering the 
co-addition. 



4.3. Photometry 

Besides the astrometric alignment, special care has to be 
taken in the photometric calibration. First, we usually 
want to co-add data obtained under different photometric 
conditions (for instance if a target is observed during dif- 
ferent nights of an observing run). Second, the individual 
images can have intrinsically different zero points^^. It is 
possible to arrive at an internally consistent photometric 
system for all exposures by comparing instrumental mag- 
nitudes of overlap sources. The following procedure of a 
relative p hotometric calibration is very similar to that de- 
scribed in iKoranvi et al.l l)l998l) . It is implemented in the 
LDAC relphotom programme. Having a full astrometric 
solution and hence an accurate table of paired objects of 
all SCIENCE images at hand, a relative photometric so- 
lution is straightforward. Given two overlapping images k 



The sky-background equalisation described in Sect. 13.61 
should remove zero point differences within an exposure but 
we do not rely on this assumption here. 
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Fig. 12. Shown are measurements for qf^ from 152 stars lying 
within one chip of 30 WFI@2.2m exposures. The quantities are 
plotted in arbitrary units. The co-addition was performed with 
SWarp and the resampled images (according to the astrometric 
solution from Astrometrix) have been used for the calculation 
of (see Sect.inj. All stars have at least 10 contiguous pix- 

els with lOo" above the sky-background in the individual images 
and hence all the measurements involved have high S/N. The 
left panels show the result for qfj^, where the co-added mosaic 
was stacked according to the astrometric solution provided by 
Astrometrix. We checked that the outliers showing a higher 
value in the co-added image are uniformly distributed over the 
images. For the right panels we shifted all individual images 
with a random offset within an image pixel before co-adding 
them. This mimics a co-addition by a simple integer pixel shift 
and we see a clear, systematic offset for all three components. 

and j, consider all i ~ I...N common objects and calculate 
the mean deviation of magnitudes K and J 



M, 



(2) 



with Wi = (c^ + where the a are the measurement 

errors of the corresponding magnitude estimates. Objects 
deviating in Ki — Ji more than a user-defined threshold are 
excluded from the following statistics. The relative zero 
points ZPi for all N overlapping images are estimated by 
a minimisation with respect to ZP^, 



X 



N 

E 

k.j 



[Mk,j - (ZPk - ZPj)Y 



(3) 



In addition we demand that J^k = to make the sys- 
tem non-degenerate. We test the robustness of our photo- 
metric algorithms with simulations and on real data in the 
next sections. We note however that our relative photo- 
metric calibration depends on two important assumptions 
that may not be fulfilled: 



Fig. 13. The same as Fig. [121 for X'''^ in arbitrary units. The 
agreement between the measurement in the mosaic and the 
expectation is not as good as for q and X™ (see Fig. 114^ and 
a systematic offset of about 1% is seen in the 11 component. 
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Fig. 14. The same as Fig. ^for X''™ in arbitrary units. For 
this quantity, the agreement between the values in the mosaic 
and the predictions is excellent for all three components. 



We assume that the relative zero points are constant 
on the scale of each image. As reported bv lKoch et alJ 
(|2004), zero point variations of up to « 0.06 mag 
within a WFI@2.2m image can occur due to non- 
uniform illuminations (e.g. scatt ered light or reflec- 
tions inside the instrument). iManfroid et all (|20^ 
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Fig. 15. Shown is the image seeing (half-light radius) of a 
set of 30 WFI<a!2.2m exposures. The mean is indicated by the 
dotted line. The solid line shows the image seeing in the co- 
added image, using Astrometrix for the astrometric solution 
and SWarp for the co-addition. Assuming that the stellar light 
profile in the individual exposures is a Gaussian, one would 
expect an image seeing as the one indicated by the dashed 
line, provided that both the astrometric solution and the co- 
addition method work perfectly well. One can see that we are 
very close to the optimal value, which is only 0.026 « l/40th of 
a pixel smaller than the actual measured image seeing in the co- 
added mosaic. One pixel corresponds to 0'.'238. This analysis 
is based on 29 stars that were in common in an arbitrarily 
selected CCD for all 30 WFI@2.2m exposures. 



show that these deviations are not corrected for by 
standard flat-field techniques (such as the apphca- 
tion of SKYFLATs or SUPERFLATs). Indeed, they 
argue that the aim of a flat sky-background is in- 
compatible with the aim of homogeneous photometry 
on the scale of an image. Similar effects have been 
observed with the CFH12K @CFHT instrument (see 
iMagnier fc Cuiflandrell2004 '). In other words, the main 
reason for the observed variations of the photometric 
zero point is the mixing of multiplicative (flat-field) 
and additive effects (scattered light, reflections) in the 
flat-fields. Simple division by such a flat then leads to 
a seemingly flat sky-background, but introduces the 
observed variations in the photometric zero point. A 
solution would be to disentangle the multiplicative and 
additive component. At the moment it is not clear how 
this could be achieved in the general case. 
We neglect that objects of different spectral type can 
be affected differently by varying photometric condi- 
tions (e.g. absorption by clouds). At this stage we have 
not investigated this effect. 



For the moment we do not take into account these effects 
in our photometric calibration. 

For a full photometric calibration of our SCIENCE 
images we need to include the results of our standard star 
calibration (Sect. 13.11(1 in the estimation of our relative 
zero points. We are currently implementing and testing 
this. 

4.4. Simulations for investigating ttie robustness of tine 
photometric calibration 

The simulations described in this section aim at quantify- 
ing the following issues of our photometric calibration: 

— How accurately can we recover relative photometric 
zero points between individual images? 

— We subtract the sky-background from our images be- 
fore they are flnally co-added (see Sect. El for details 
on the procedure). Does this procedure introduce any 
systematic bias in the flux measurements of faint ob- 
jects? 

To this end we simulate d WFI exposures w ith Emmanuel 
Bertin's skymaker^^ (see lErben et al.l200ll) and processed 
them through our pipeline in exactly the same way as real 
data. In this way we include all possible uncertainties from 
measurements, possible pixelisation effects etc. in our sim- 
ulations. The simulations have the following characteris- 
tics: 

— We created 25 exposures of 8 chips in the layout of the 
WFI@2.2m telescope (i.e. a pixel size of 0'.'238 and a 
chip size of 2k x 4k) . Each exposure mimics a 500 s in- 
tegration in an R broad-band filter. The images do not 
contain any flat-field or fringe residuals but resemble 
perfectly pre-reduced data. 

— Each exposure is put randomly within a 2' x 2' dither 
box. No geometric distortions have been introduced. 

— Each exposure has a random image seeing between 0'.'9 
and I'.'l (no PSF anisotropy) and a random sky i?-band 
background between 20.0 and 20.5 mag/arcsec^. 

— Each image is assigned a random magnitude zero point 
between 24.9 ± 0.3, i.e. there is no correlation of zero 
point with exposure, extinction or gain. 

The images were processed through our pipeline starting 
at the astrometric calibration up to the flnal mosaic. In 
addition we simulated with skymaker a single 10k x 10k 
exposure with a formal integration time of 25 x 500 s (sky- 
maker mosaic in the following) for flux comparisons with 
the co-added mosaic. The accuracy of the zero point esti- 
mation is shown in Fig. 1161 We see that our procedure re- 
covers the input values with a formal uncertainty of 0.002 
mag. In Fig. 1171 we show the comparison of flux measure- 
ments from the co-added mosaic with the skymaker mo- 
saic. No systematic deviations in the magnitudes up to 



the TERAPIX tool for the simulation of astronomical im- 
ages 
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m.0C| mag 



Fig. 17. Investigation of possible fiux measurement biases in our mosaics: We simulated 25 dithered WFI exposures with 
skymaker and processed them through our pipeline in exactly the same manner as real data. In addition we created one large 
skymaker mosaic not undergoing any pipeline processing. The figure shows the difference of magnitude estimates in these two 
images, i.e. it mainly reveals possible biases due to inaccurate sky subtraction in the mosaic. The left panel shows the raw 
differences, the right panel bins with a width of 1 magnitude. We notice that our processing does not introduce systematic 
biases up to the faintest magnitudes. See the text for more details. 



the faintest objects are visible. We conclude that our pro- 
cedures accurately recover relative zero point deviations 
between images and that our mosaics do not suffer from 
systematics in estimating magnitudes. Errors connected 
to the algorithms are, at least for WFI@2.2m, small com- 
pared to systematic errors (variations of about 0.05 mag 
within an image) due to possible CCD illumination prob- 
lems (see above). 



four fields with the Warsaw 1.3 m telescope at the Las 
Campanas Observatory. 

Using the data we solved the following transformation 
where we neglect the airmass-term which is constant over 
the field-of-view: 

b = B + bi-{B-V) + b2 (4) 
V = V + vi-{B-V)+V2. (5) 



4.5. Photometric measurements in NGC 300 

To finally quantify the quality of our current photometric 
calibration, w e compare flux m easurements in the field of 
NGC 300 fsee lSchirmer et al.l[2003.) with previous l y pub - 
lished secondary standards bv IPietrzvhski et"al] |2002|) . 
PSF photometry was carried out on co-added B (expo- 
sure time: 39 600 sec) and V (exposure time: 37440 sec) 
WFI @2.2m images of NGC 300 by using the DAOPHOT 
task llStetsonlllOS ?) implemented in IRAF Version 2.12.2. 
We used about 100 out of 390 secon dary photometric 
standa rd stars that were published by IPietrzvhski et al.l 
l)2002|l to uniformly cover our field-of-view. Areas, in 
which crowding became significant or where estimations 
of the sky-background became difficult (galaxy centre, spi- 
ral arms, saturated stars), were excluded. Fig. 1181 shows 
the spatial distribution of the chosen standard stars in 
the co-added image. The galaxy-centre of NGC 300 co- 
incides with the field centre. T h e pho tometric standards 
published in IPietrzvhski et all l)2002fl were observed in 



Capital letters represent standard star magnitudes in 
B and V respectively whereas lower case letters indicate 
the measured instrumental magnitudes. Values of &[i,2] 
and V[i 2] have been determined by solving these equations 
for the standard stars using the Levenberg-Marquardt al- 
gorithm implemented in the IRAF routine INLFIT. This 
leads to values of bi = -0.28 ± 0.007, 62 = 0.69 ± 0.006, 
vi = 0.078 ± 0.005 and V2 = 0.89 ± 0.005. In Fig. Hi we 
plot the residuals of this fit as a function of _B-band (left- 
hand side) and l/-band (right-hand side) magnitude and 
see that there is (to a first approximation) no significant 
tendency for the residuals to correlate with magnitude^^. 
Calculating the standard deviation a for these fits we get 
(t{B) — 0.037 and a{V) — 0.029 with maximum residuals 
of about 0.07 magnitudes. 

This means the errors and/or uncertainties introduced 
by a global solution for the majority of stars are smaller 
than 0.04 magnitudes and are therefore negligible for most 
photometric studies taking into account that they are of 



this especially holds for stars fainter than mag=19 
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Fig. 19. Residuals vs. magnitude plots of 107 chosen standard stars. The left panel shows the results for the i?-band, 
whereas the right panel shows those for the V^-band 
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Fig. 16. The figure shows the error distribution of relative 
photometric zero point recovery in Monte Carlo simulations. 
The a;-axis displays the difference Am (in magnitudes) between 
input and recovered relative zero point, the y-axis the recovery 
probability. Our algorithm recovers the input values formally 
with Am = ± 0.002. See the text for more details on the 
simulations. 



Fig. 18. Position of the chosen standard stars in the co- 
added y-band image. Crosses denote the positions of the 
secondary standards from Pietrzvhski et al. (2 00^ . the 
big dot the centre of NGC 300 and the f illed h exagons 
two photoelectric standards from lGrahamI (^^). 



the same order as the uncertainties for some of the sec- 
ondary photometric standards. 

We also tested the global solution in relation to the 
different chips. For this, we fitted the global solution sep- 



arately to areas where mostly images of a particular chip 
overlap and checked the standard deviation of these fits. 
The results are presented in Table 0] and it becomes ob- 
vious that the standard deviation for a single chip differs 
from the standard deviation of the global solution only 
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Table 4. Standard deviation of the global solution ap- 
plied to each chip. The relation between chip number and 
actual position on the mosaic can be seen in Fig. ^ 



Chip No. 


a(B) 




1 


0.035 


0.024 


2 


0.040 


0.034 


3 


0.039 


0.033 


4 


0.040 


0.031 


5 


0.030 


0.026 


6 


0.036 


0.033 


7 


0.039 


0.033 


8 


0.039 


0.027 


1-8 


0.037 


0.028 



by a constant smaller than 0.005. Due to this, our global 
solution can be regarded as sufficient and therefore shows 
that our photometric calibration is capable of removing 
chip-to-chip variations in a more than satisfying manner. 
This has also to be seen in the context that the observa- 
tions of the very extended NGC 300 galaxy were spread 
over 34 nights in six months, encompassing very different 
sky conditions. 

Finally we have performed a check of our transfor- 
mation with_two photoelectric standards published by 
which are not saturated in our images. 
Applying the global solution to these photoelectric stan- 
dards leads to residuals of —0.036 and —0.069 in V and 
—0.055 and —0.022 in B. Considering the fact that on the 
one hand the photometric errors for these stars are about 
0.05 magnitudes, and on the other hand some of our fit- 
ted standards had residuals up to 0.07 magnitudes, these 
residuals are in good agreement with the previous results. 

We point out that the checks performed so far give 
us confidence that to date the relative photometry of the 
pipeline permits us to reduce and to co-add images to a 
photometric accuracy of about 0.05 magnitudes (the av- 
erage between the standard deviations of the fits t o sec- 
ondary photometric standards of Pietrzynski et al. ,120021 
and photoelectric standard stars from Graham, Il98ll) . 
This is in agreement with the expected errors due to the 
above mentioned illumination issues. 



4.6. Set quality characteristics 

The object catalogues generated for the astrometric and 
photometric calibration give us a good overview of the 
data quality at this stage. They allow us to investigate 
in more detail the PSF properties of individual exposures 
which is crucial in weak lensing studies. Also the night sky 
conditions (transparency, sky brightness) can be studied 
in detail and images that should be excluded from the 
final co-addition can be identified. Fig. [50] discusses in de- 
tail PSF properties of WFI@2.2m and Fig. ED shows an 
overview of our currently implemented quality control. 



5. Weighting and Flagging 

With astrometric and photometric information for all im- 
ages at hand we can finally co-add them and produce a 
deep mosaic. This section motivates our choices for the 
implemented co-addition procedure which i s very close to 
that o f the EIS Wide Survey described in iNonino et alJ 

As long-exposed SCIENCE frames are dominated 
by Gaussian sky noise the optimal result (in terms of 
noise properties in the co-added image) is obtained by a 
weighted mean co-addition of all input images. However, 
this is not straightforward as many image pixels carry non- 
Gaussian noise properties (such as bad columns, vignetted 
image regions, cosmics) or other defects that we would like 
to exclude from the co-addition process (such as satel- 
lite tracks or extended stellar reflections). To perform a 
weighted mean co-addition all bad image pixels have to 
be known beforehand and assigned a zero weight in the 
co-addition process. We note that identifying defects be- 
fore the co-addition has several advantages over rejecting 
them during the stacking process: 

— Rejection methods during the co-addition, such as me- 
dian stacking or sigma-clip rejection require a mini- 
mum number of 5-10 input images for robust statistics. 
With prior knowledge of the defects good results can 
already be obtained with a small number of images. 

— Co-additions based on median or sigma-clip rejection 
are problematic when observations under varying see- 
ing conditions have to be combined. In this case object 
profiles in the co-added image might be affected by the 
rejection algorithm. 

— A practical advantage is that a weighted mean with- 
out rejections is a strictly linear process. Hence, new 
images can directly be added to an already existing 
mosaic and not all input images need to be co-added 
again. 

We identify bad pixels in the individual exposures as fol- 
lows: 

1. Every chip contains hot and cold pixels, i. e. pixels 
that always have a high/low charge value. They are 
typically stable over the period of an observing run 
and are most effectively identified in the masDARK 
frames (see Fig. l22|l . If DARK frames are not at hand, 
SKYFLATs and/or SUPERFLATs often provide a suf- 
ficient bad pixel map. 

2. Saturated image pixels are identified by applying a 
pixel value threshold to the SCIENCE frames. 

3. Pixels affected by cosmic rays are not identified as eas- 
ily as their appearance is different in each SCIENCE 
frame. We detect them with SExtractor with a spe- 
cial filter for cosmic ray detection generated with EYE 
l)Bertinll200lt iNonino et al. 1999 ). The filter we actu- 
ally use was specifically designed for cosmic ray re- 
moval on EMMI@NTT data in the framework of the 
EIS Wide project. It is difficult to say for which other 
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Fig. 20. PSF properties of WFI@2.2m fields: For each image within a set we characterise PSF properties (PSF size and 
anisotropy) over the complete field-of-view of a WFI. To this end we run the KSB algorithm [a weak lensing analysis technique 
for accurate shape measurements of astronomical sources; see lKaiser et al.l l^^^)] over selected stars of all images. The upper 
left panel shows PSF anisotropies for an int ra-focal , the upper right for a focal and the lower left for an extra-focal exposure. 
The sticks represent amplitude (given as \/ef-l-ef) and orientation (f> — 0.5 arctan(e2/el) of the PSF anisotropy. The chosen 
scale for the stick length is the same for these three plots in order to show the increase in the anisotropies with respect to the 
focused exposure. The mean stellar ellipticities are 0.066, 0.09 and 0.059, respectively. The lower right panel depicts typical 
PSF anisotropies of a WFI@2.2m i?-band mosaic ('^ 57 exposures with ~ 500 s exposure time each). Note that the largest PSF 
anisotropies in the mosaic are as small as ~ 0.01. Compared to the other three PSF plots, a different scale for the stick length 
was used in order to clearly show the anisotropies. We note that WFI@2.2m has excellent PSF properties if data are obtained 
under favourable focus conditions. 



instruments it can be used without modification or to 
quantify its actual performance. The visual impression 
from our co-added WFI@2.2m frames is that the de- 
tection efficiency for small-scale cosmic ray features is 



very high and gives a sufficient image cleaning for this 
camera (see Fig. 
4. Other extended defects such as satellite tracks or re- 
flections from bright stars are masked by hand. They 
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Fig. 21. Shown are quality control plots for a set of 28 R band exposures from WFI@2.2ni. They primarily allow us 
to select the images that enter the final co-addition process. The two upper left plots show the seeing distribution 
in the set, where the seeing value of an exposure is estimated from the mean seeing in the individual chips. The 
distribution of seeing variations within an exposure is also shown (second left plot). The two lower left plots show 
the night sky amplitude (normalised to an exposure time of one second) during the observations. They allow us to 
identify exposures taken during twilight or unfavourable moon phases. In the upper right we have a quick check on 
the photometric conditions through the distribution of relative zero points (see Sect. 14. Hi . The standard deviations of 
relative zero points show that our sky-background equalisation (see Sect. l3.6|l adjusts zero points typically to 0.01-0.03 
mag within a WFI@2.2m mosaic. In the lower right we see the PSF ellipticity distribution (see Fig. I20() . Here each 
point represents the mean ellipticity value from all stars within an exposure. 



have been identified in the manual pass through the 
data during the Run processing (see Sect. 13.1011 . 

In this way we can generate a FLAG map for each im- 
age. It is an integer FITS image containing zeros for every 
good pixel and values characterising the different defects 
otherwise. We use these FLAG maps in the creation of 
WEIGHT maps (see below) and to mark standard stars 
whose magnitude measurement is affected by a defect (see 
Sect. EH. 19 

Besides to the actual image co-addition algorithm we 
need to pay special attention to the noise properties in the 
final co-added image. The gaps in multi-chip cameras and 
the different intrinsic gains of the various detectors lead to 
complex noise properties in a co-added image consisting 
of dithered exposures. These noise variations need to be 
taken into account properly when estimating a threshold 
for object detection or for estimating errors based on pixel 
noise (such as object fluxes for instance). For a full char- 
acterisation of the noise, another image besides the final 

The FLAG and WEIGHT maps described in this section 
are actually created on the Run level before the standard field 
processing. Because of their very close relation to the image 
co-addition we shifted their description to the Set processing. 



co-added image, the WEIGHT map, describing the rela- 
tive noise properties for each image pixel, is created during 
the co-addition process. We arrive at this WEIGHT map 
in the following way: 



1. We produce a WEIGHT map for every input image 
to be co-added. The starting point for these individ- 
ual WEIGHTS is the masFLAT that is rescaled to a 
mode of unity. It provides information on the relative 
sensitivity (and hence noise) variations within an im- 
age. This basic map is modified by the FLAG map by 
setting the value of defect pixels to zero. 

2. These individual WEIGHTS are co- added together 
with the science images. The co-added WEIGHT pro- 
vides a full characterisation of the relative noise prop- 
erties later. The exact co-addition procedure is de- 
scribed in Sect.lHl Fig. 1241 shows an example for a co- 
added WEIGHT map and Fig.|^its merits for object 
detection. 



FLAG and WEIGHT maps are produced with the 
TERAPIX tool weightwatcher fsee lBertia.lQQS,) . 





Fig. 23. The left panel shows 25 co-added WFI@2.2m frames where no cosmic ray cleaning has been performed. The 
right panel shows the result with the cleaning procedure turned on. We find that for WFI@2.2m frames SExtractor 
(in connection with a neural network filter generated with EYE) performs a very efficient detection of small-scale cosmic 
ray features. 




Fig. 24. Shown is the WEIGHT map (left panel) characterising the noise properties of a co-added science image and 
a zoom- in to its centre (right panel). The lighter the colour, the higher the relative weight (the lower the noise) of the 
pixel. The darker regions at the positions where different chips meet have about half the weight of well-covered regions. 
Different weights between regions where the same number of input images have contributed show intrinsic sensitivity 
variations. We note that the noise structure is quite complex and cannot be taken into account appropriately without 
the WEIGHT map. These maps are used by SExtractor in the object detection phase and for noise calculations. 
Having the map at hand one also does not need to cut the outermost regions of the coadded images where the noise 
is considerably higher than in the inner regions. See also Fig. 1251 



6. Image co-addition 

Before the final co-addition, all images are brought to the 
same background level by subtracting the night sky. To 
accurately model the sky-background and to avoid biased 



estimates close to large astronomical sources we apply the 
following two-stage process: 

1. We run SExtractor and detect all large-scale objects 
having at least 50 adjacent pixels with 1.5ct over the 
sky-background (these values are for WFI@2.2m and 
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Fig. 25. SExtractor source catalogue obtained from the same mosaic (with an area of 95" x 95" ) without (left) and with 
(right) a WEIGHT image as an additional SExtractor argument. The number of spurious detections in regions with higher 
noise is obvious when no WEIGHT image is used as an additional input. 




Fig. 22. Shown is a masDARK image and a zoomed in 
region. In a stack of several long-exposed DARK frames 
hot and cold pixels, which often come in groups and affect 
complete rows/columns, show up with high significance 
and are easily identified by applying a pixel value thresh- 
old to the masDARK. 

can vary according to the data set). All image pixels 
belonging to detected objects are replaced with the 
image mode. 

2. From this object-subtracted image we create 
a SExtractor BACKGROUND check-image 



(BACK_SIZE=90 for WFI@2.2m images) and 
subtract it from the original SCIENCE image. 

We showed in Sect. l4.4l that this sky subtraction does not 
introduce a significant bias in the object photometry on 
the CO- added images. 

For the co-addition we can use Richard Hook's EIS 
Drizzle which is implemented within IRAF or SWarp 
(E. Bertin, TERAPIX). Both software packages perform 
a weighted mean co-addition on a sub-pixel level taking 
into account the full astrometric solution and allowing 
a variable output pixel size. E IS Drizzle was dev eloped 
within the EIS Wide project i Nonino et aD llQQO'l to co- 
add quickly large volumes of data with low demands on 
CPU and virtual memory. It uses strongly simplified cal- 
culations during the pixel resampling process and allows 
only strictly linear co- additions (see below). In contrast, 
SWarp offers a large variety of pixel resampling and co- 
addition algorithms. 

SWarp first undistorts and resamples all input 
SCIENCE and WEIGHT images according to the astro- 
metric solution. The user can choose between several so- 
phisticated kernels for the pixel re mapping. We use by 
defauh the LANCZ0S3 kernel (see iBertinI I2OO2I for de- 
tails on SWarp ' s resampling kernels). The final co-addition 
of the resampled images is performed in a second pass. 
Having all resampled input images belonging to a given 
output pixel at hand simultaneously, SWarp can calcu- 
late the final result in a variety of ways, such as me- 
dian, mean or the weighted mean which is our method 
of choice. EIS Drizzle's implementation is significantly 
different. It uses a method known as forward mapping. It 
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performs the complete co-addition in one pass by putting 
the input images consecutively on the output grid. Hence, 
EIS Drizzle is limited to linear co-addition methods. 
Applying the astrometric solution, an individual input 
pixel is mapped somewhere on top of several neighbour- 
ing output pixels, and is distorted and rotated. Its flux 
is then distributed accordingly amongst these output pix- 
els. The EIS Drizzle approach strongly simplifies the 
calculation of the flux distribution, in the sense that only 
non-integer shifts are taken into account, whereas rota- 
tions and distortions of the mapped pixels are neglected. 
This corresponds to the Turb o Kernel in the curren t 
Mult iDr izzle implementation llKoekemoer et alJ l2002|) . 
See also iFruchter &: Hood 1)2002!) for a more detailed de- 
scription of the drizzling approach. We discuss the differ- 
ences of the pixel resampling kernels of SWarp and EIS 
Drizzle in more detail at the end of this section. 

Choosing a weighted mean co-addition for SWarp, both 
software packages calculate the value of an output pixel in 
the same way. Four factors contribute to the final result 
in our weighted mean co-addition. Given are the value 
li, representing the part of an input pixel that goes to 
a specified output pixel /out in the co-added image. Wi 
represents the associated value in the WEIGHT map. li 
is scaled with factors fi to the consistent photometric zero 
point and normalised to a fixed exposure time of 1 s. This 
scaling reads 



/. = 10 



-OAZPi 



(6) 



where ti is the exposure time and ZPi the relative pho- 
tometric zero point. All images are in addition weighted 
according to their sky noise. This weight scale is given by 



f2 

^ sky jiJ i 

We take into account that the noise also scales with the 
flux scale fi. The values /out and Wout in a stack of N 
exposures then read 



laut — 



N 

Wont = ^'^^ 
i=l 



(8) 



Besides the co-added SCIENCE and WEIGHT mosaics 
EIS Drizzle can produce a CONTEXT map which allows 
the identification of all input frames that contribute to a 
given pixel in the co-added images. 

As standard sky projection of the mosaic we use the 
TAN projection fsee ICreisen fc Calabrettall2002l for fur- 
ther information on sky projections) with an orientation 
such that North is up and East to the left. A reference 
coordinate can be specified for the co-addition, associated 
with an integer reference pixel. Thus, if multi-colour infor- 
mation is available for a particular Set, the mosaics in the 
different bands can be registered with sub-pixel accuracy 
if required. 

The user can vary the fraction of the input pixel entering 
this procedure (drizzle PIXFRAC parameter). We always use 
the original pixel size in our co-additions (PIXFRAC=1). 



In the following we show a comparison of photomet- 
ric properties from objects extracted from the same set 
co-added with EIS Drizzle and SWarp. We expect that 
the simpler resampling kernel from EIS Drizzle leads to 
stronger noise correlations and image blurring than the 
more sophisticated SWarp approach. This is illustrated in 
Fig. HSl and Fig. |13 The EIS Drizzle co-addition has 
a slightly larger image seeing and SExtractor chooses a 
larger optimal radius for its MAG_AUTO parameter. In 
contrast, EIS Drizzle has no significant impact on the 
PSF anisotropics (see Fig. Also more sophisticated 
weak lensing analyses (such as cluster mass reconstruc- 
tions) revealed that results based on object shape mea- 
surements do not differ significantly between the two co- 
additions. 
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Fig. 26. Comparison of the image seeing between two mo- 
saics created with SWarp and EIS Drizzle. Shown are the half- 
light radii for unsaturated stars in the two co-additions. The 
swarped image has an image seeing that is 0.09 pixels smaller 
than the one for the drizzled image. The mean values for the 
image seeing are 0'.'95 for the drizzled image and O'.'Ql for the 
swarped one. The same astrometric solution was used for both 
co-additions. Thus EIS Drizzle slightly increases the size of 
the PSF, an effect of its simplified kernel. In Fig. l28l it is shown 
that the PSF anisotropics are identical for both co-addition 
strategies. 



7. First Data quality assessments of co-added Sets 

The implementation of a thorough, automatic quality 
check on our co-added mosaics is still in its infancy. In 
this section we describe basic tests on extracted object 
catalogues that are done without any user interaction. 
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Fig. 27. We show comparisons of magnitude estimates from an EIS Drizzle and a SWarp co-addition from the same set. 
Both measurements are in very good agreement if the flux is measured within a fixed aperture (here with a diameter of 3'.'5). 
However, for SExtractor's MAG_AUTO estimate (Kron like total magnitudes), the aperture used for faint sources is larger in 
the drizzled images than in the SWarp mosaic. The reason is the more correlated noise in the EIS Drizzle approach. 



7.1. Galaxy counts 

The programme SExtractor is used to create a raw cat- 
alogue of all objects that consist of at least 3 contiguous 
pixels (DETECT_MINAREA=3) with a flux 2a above the 
flux of the sky-background (DETECT_THRESH=2). The 
source extraction is done on a filtered image; we use a 
normalised Gaussian filter with a full width half maxi- 
mum of 4.0 pixels (FILTER_NAME=gauss_4.0_7x7.conv). 
This conservative threshold is chosen in order to min- 
imise the number of spurious detections. In the follow- 
ing, we use the SExtractor parameters MAG_AUTO for 
magnitudes and FLUX_RADIUS for the half-light radius. 
To create a galaxy catalogue, the SExtractor parameter 
CLASS_STAR in combination with the half-light radius is 
used; saturated objects are rejected. We define every ob- 
ject which has CLASS_STAR less than 0.95 as a galaxy. 
A simple check of this selection is a magnitude over half- 
light radius plot. All stars have the same half-light radius 
and therefore show up as a vertical branch in this plot, 
see Fig. EHl 

We count the number of galaxies in 0.5 mag wide bins per 
one square degree. To normalise the area to one square 
degree, we take into account that each object occupies an 
area in which fainter objects can not be detected. For this 
correction we use the SExtractor parameter ISOO, which 
is the isophotal area above the analysis threshold. We note 
that in the case of empty fields this effect is almost negligi- 
ble. An error-weighted linear regression to the logarithmic 
galaxy counts is performed and the slope, dlog N/ dmag 
determi ned. We routi nely compare our galaxy counts with 
those of iMcCrackeiTet al.. (,2o6i) . see Fig. EOl With this 
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Fig. 29. Magnitude vs. half-light radius plot. Grey 
crosses: stars; appear as a vertical branch in the plot and 
are selected by the CLASS_STAR parameter and half-light 
radius. Dots: extended sources. 



comparison, a rough test of the magnitude zero point and 
the limiting magnitude can be performed. 



30 



T. Erben et al.: Methods for WFI data reduction 



— 1 — I — I — I — I — I — I — I — I — I — I — I — I — 1— 



t X 



-a.D4 - 



■ if *- t* 

.7- i .■^i^' T-: 



.Hi 




K ,: i - -^i-.V-t.. 1 • 



I I 



^,02 



■ 1 driula 



a.Di 



DUUft -(>.CG 




_l L 



I I 



«1 Akdip 




Cl 3010 .3CKID £l I£C1D OdOD ^MD SOHJ 

^Ei^ : li^i &fi11 ■'^ ^EiJ: assail fefilf ' ^ 

Fig. 28. PSF anisotropics in a drizzled (left column) and a swarped (right column) mosaic of the same data set. The patterns 
are virtually identical. Yet EIS Drizzle marginally increases the size of the PSF, as was shown in the left panel of Fig. 1261 The 
mean PSF anisotropy for this particular mosaic amounts to mere 0.022 ± 0.015. 



7.2. Clustering of extended sources 

A further test for the quality of the co-added image is the 
clustering of sources. For that purpose we use the two- 
point angular correlation function, lu{9), where uj{6)S0 is 
the excess probability of finding a pair separated by an 
angle between 6 and 9 + 69. We estimate this quantity by 
creating a large number of random catalogues (by default 
40 mock catalogues are created) and count the pairs within 
the data catalogue, DD, within the random catalogue, RR 



and between the data and rand om catalogues, DR. Th e 
estimator for uji9), proposed by iLandv fc Szalavl l)l993j) . 



iui9) = 



DD - 2DR + RR 
RR ■ 



(9) 



The random fields must have the same geometry as the 
data field. Therefore we calculate an obscuration mask out 
of the number density of extracted sources as follows. A 
mesh with 512 x 512 mesh cells is placed on top of a data 
field and the number of objects in each cell is counted. 
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Fig. 30. The figure displays the logarithmic galaxy counts 
in 0.5 mag bins per one square degree from our re- 
duction of the Chan dra Deep Field South (CDFS; see 
iGiavalisco et al.ll2004|) . The error bars are due to Poisson 
noise; the line fit is an error- weighted linear regression in a 
magnitude range between the saturation and the limiting 
magnitude (here: R G [16.3; 23.8]). For the normalisation 
of the area we take into account that each detected object 
occupies an area in which fainter objects cannot be de- 
tected. As a comparison to our galaxy counts we also plot 
the galax y number counts from the CFH12K-VIRM0S 
deep field ((McCracken et al.ll20n3l). For the CDFS we also 
plot the number counts from lArnouts et al 

Then the cell count matrix is smoothed with a Gaussian 
kernel of about 512/60 cells FWHM, and all matrix ele- 
ments lower than a given fraction (we use 75% as default 
but varying this parameter between 50% and 80% does 
not change the results significantly) of the mean number 
of galaxies inside a cell are defined as a masked region. 
The borders of a field, bright stars, and part of their halos 
are masked automatically by this method (see Fig. I31|l . 

To maximise computational speed, we perform our cal- 
culations by creating an index tre e for galaxy position as 
explained in IZhang fc Pe"^ l)2004|) . The error bars in our 
check-plots for each angular bin are simply estimated by 
Poisson noise and are therefore a lower limit to the uncer- 
tainty in Lu{9). 

As an example we present the two-point angular cor- 
relation function of galaxies in our i?-band reductions for 
11 WFI-fields in the EIS Deep Pubhc Survey (EIS DPS). 
For three different magnitude bins we fitted a power law, 
tij{9) = A9^ — C, to the data. The variable C is the so- 
called integral constraint (see for instance iRoche et alJ 
which only becomes important for large fields at 
larger scales {6 > 2'). We therefore neglect the integral 
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Fig. 32. The two-point angular correlation function of 
galaxies for different magnitude bins, for 11 WFI-fields 
in the DPS. We show the logarithm of the amplitude, 
log uj{0), as a function of the angular separation 9 in ar- 
cmin. Here the error bar in each bin is due to the field-to- 
field variance. We perform a simple error-weighted linear 
regression in the angular interval 9 e [0'.'06; 2' ] (where the 
integral constraint is negligible) to determine the slope 5 
and the amplitude Ai = uj{9 — 1'). 



constraint and perform the fit for small scales {9 < 2'). 
The results are shown in Fig. 

The results from our clustering analysis can be cross- 
checked by considering the aperture number count disper- 
sion {N^){9). It is directly related to the angular correla- 
tion function to (9) by 



(10) 



where the function is defined as eq. (35) in 

ISchneider et al.l l)2002() and it has the nice property of be- 
ing independent of the integral constraint. As an exam- 
ple we use the measured lo{9) for the magnitude interval 
R e [19.5; 21] from the 11 WFI-fields of the last section 
to calculate (iVobs)(^)> see Fig. ESI To compare the slope 
5 and amplitude Ai of the fitted power law to the mea- 
sured Lo{9) of the previous section (Fig. {N\^){9) is 
calculated for the power law, ijj{9) = Ai 9^ . The results 
show that the slope S obtained from the power law fit 
of the angular correlation function is correct. The ampli- 
tude Al, however, seems, with Ai — 0.133, to be a bit 
larger; this small discrepancy can easily arise due to the 
fact that {N^hs)i^) is calculated using lu{9) over the en- 
tire 6'-range and that the w-fit is slightly influenced by 
the integral constraint. As a comparison. Fig. displays 
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Fig. 31. Obscuration mask for the CDFS. Left panel: Distribution of detected sources. Right panel: FITS image of 
the CDFS. White contours indicate masked region. In deep fields, our automatic obscuration mask based on object 
density variations reliably marks large-scale astronomical sources and noisy borders which would significantly infiuence 
the area of our correlation function analysis. This approach turned out to be sufficient for quality check purposes. 



the function (N^) (6) for the fit parameters determined by 
[llockstra et al. ( 2002) for i?c-band data of the same mag- 
nitude range. Note that this can only be a rough compari- 
son, because of different source extraction algorithms and 
slightly different i?-band filter. 

8. Conclusions and Outlook 

We have presented our image processing methods for 
multi-chip cameras that we developed in the course of the 
GaBoDS project. A significant fraction of GaBoDS is a 
virtual survey in which observational data were collected 
from the ESQ Science archive within an ASTROVIRTEL 
program^^ primarily for weak gravitational lensing stud- 
ies. This allowed us to test and to apply our procedures 
on data sets which were acquired for a large variety of 
scientific programmes (e.g. deep field observations, search 
for moving objects, Cepheid studies in nearby galaxies) 
and obtained in many different ways (very compact to 
very wide dither patterns. Sets observed within a single 
night or over several years). Our experiences regarding 
the processing of these data sets from WFI@2.2m are the 
following: 

— The techniques described perform very well on empty 
field observations with WFI@2.2m. The excellent op- 
tics of this instrument allow an accurate astrometric 
alignment of images which is cru c ial fo r weak lens- 
ing studies l|Schirmer et all 120031 12004 - Also very 



ASTROVIRTEL Cycle 2: Gravitational lensing studies in 
randomly distributed, high galactic latitude fields; P.L Erben 



0.03 



0.02 



0.01 



■A, = 0.133; (5 = -0.9 
-A, = 0.121; (5 = -0.9 
A, = 0.115; (5 = -0.7 
(Hoekslra et al. 02) 




10 

6 [arcmin] 



15 



20 



Fig. 33. The aperture number count dispersion (N^bg){0) 
as a function of angular scale 6 in arcmin for 11 WFI- fields 
in the magnitude interval R e [19.5; 21]. The error bars on 
the function {N'^){6) are due to the field-to-field variance 
of the 11 fields. Note that the points are correlated. The 
lines display the function (N^^) (9) assuming a power law 
for the angular two-point correlation function of the form 



w{9) = Al 9 for different parameters Ai 



1') 



and 6. Solid line: {N'^^) (9) calculated for the fit parameter 
obtained from the angular correlation function (see Fig. 
I32|l : dotted line: comparison with (N^)(9 ) calculated for 
the fit parameters determined bv iHoekstra et al.. (,200^ : 
dashed line: best fit to (N^){9). 
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large data sets, obtained over several years, can be 
processed. For example, we collected and reduced all 
WFI@2.2m observations from the CDFS, which con- 
sist of more than 100 individ ual exposures in each 
of the U, B, V and R bands jGiavalisco et all 120041: 
iMobasher et al.ll20o4 . 

— Data sets from crowded fields or from large-scale ob- 
jects, whose extent is comparable to the field-of-view, 
can be processed with good res ults but require sub- 
stantial manual intervention fsee lSchirmer et aDl200,'4 
for processing details of the field around NGC 300) . 

— Our absolute photometric calibration is currently ac- 
curate to about 0.05 mag as discussed in Sect. 14.31 
and Sect. 14.51 Also the comparison of our standard 
star calibrations with independent measurements in 
Sect. 13.111 shows errors of the same order. However, 
it seems that non-uniform illuminations that probably 
contribute most to the error budget, do not strongly 
depend on wav elength and hence errors on colours 
are smaller (see iKoch et al.ll2004l) . Our reduction of 
the CDFS data gives good resul ts in photometric red- 
shift studies of different groups llMobasher et alJl2004l: 
iGabasch et alJl2004 Iffildebrandt et alJl2004|) . 

As was described in Sect. |21 our processing system has 
already been used with a larger variety of single- and 
multi-chip cameras. However, more sophisticated tests on 
the astrometric and photometric accuracies, as presented 
here for WFI@2.2m, have not b een done except for the 
single-chip camera F0RS1@VLT (iBradac et al, 2004) an d 
F0RS2@VLT with a two-chip camera l)Erben et al.l20"o^ . 

We are currently performing a thorough comparison of 
released fields from the EIS DPS with our own reductions 
of this survey. This will give us additional insights in the 
quality and the properties of our reductions and allows us 
to identify useful and necessary quality tests for the co- 
added mosaics. This work will be presented in Dietrich et 
al. (in prep.). 
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Appendix A: Pipeline Image header 

Our pipeline replaces the original FITS headers of all in- 
dividual CCDs by a new one containing only a minimum 
set of keywords. In this way we unify the headers for all 
instruments and avoid inconsistencies especially in the as- 
trometric calibration. An example header for WFI@2.2m 



is: 






SIMPLE = 


T 


/ 


BITPIX = 


16 


/ 


NAXIS 


2 


/ 


NAXISl = 


2142 


/ 


NAXIS2 = 


4128 


/ 


BSCALE = 


1. 


/ 


BZERO 


32768. 


/ 



CTYPEl = 


'RA— TAN' 


/ 


CTYPE2 = 


'DEC— TAN' 


/ 


CRPIXl = 


-416. 


/ 


CRPIX2 = 


-224. 


/ 


CD1_1 


-6.61E-05 


/ 


CD2_2 


6.61E-05 


/ 


CD1_2 


0. 


/ 


CD2_1 


0. 


/ 


CRVALl = 


12.505009 


/ 


CRVAL2 = 


-52.15978 


/ 


RADECSYS= 


'FK5' 


/ 


FILTER = 


'BB#Rc/162_ESD844' 


/ 


AIRMASS = 


1.113885 


/ 


EXPTIME = 


599.9176 


/ 


EQUINOX = 


2000. 


/ 


IMAGEID = 


3 


/ 


GABODSID= 


1401 


/ 


EISID 


10 


/ 


OBJECT = 


'BPM16274_3' 


/ 


ZP 


-1.0 


/ 


COEFF 


-1.0 


/ 


DUMMYO = 





/ 


DUMMY 1 = 





/ 


DUMMY2 = 





/ 


DUMMY3 = 





/ 


DUMMY4 = 





/ 



END 

The IMAGEID, GABODSID and EISID are unique iden- 
tifiers for the chip position within the mosaic, the night 
of the observation and the image respectively. ZP and 
COEFF will finally contain magnitude zero point and ex- 
tinction coefficient. At the end we introduce 20 DUMMY 
keywords allowing the user to transfer important infor- 
mation from the original headers. For an explanation 
of the rest of the keywo rds see l|Hanisch et all I2OOII: 
iGreisen fc CalabrettalEio^. 

Appendix B: Released software 

We will release all modules of our processing system that 
are based on publicly available code. Notes on their re- 
trieval will be given in the final, accepted version of this 
paper. 



